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Abstract 



We present a review of the latest developments in one-dimensional ( ID I optical wave 
turbulence ( |OWT I . We describe theory, numerics and experimental observations for a 
|lD| optical system. This system is described by two coupled nonlinear equations, which we 
explore within two wave limits allowing for the expression of the evolution of the complex 
amplitude in a single dynamical equation. We consider a long- wave limit corresponding to 
waves with wave numbers smaller than the electrical coherence length of the liquid crystal, 
and the opposite limit, with wave numbers are larger. We show that both of these systems 



are of a dual cascade type, analogous to two-dimensional (2D) turbulence, which can be 



described by wave turbulence (WT) theory. We show that the cascades are induced 
by a six-wave resonant interaction process. |WT| predicts several stationary solutions 
(non-equilibrium and thermodynamic) to both the long- and short-wave systems, and 
we investigate the necessary conditions required for their realization. Interestingly, the 



long-wave system is close to the integrable ID nonlinear Schrodinger equation (NLSE) 



(which contains exact nonlinear soliton solutions), and as a result during the inverse 
cascade, nonlincarity of the system at low wave numbers becomes strong. Subsequently, 



due to the focusing nature of the nonlinearity, this leads to modulational instability (MI) 
of the condensate and the formation of solitons. With the aid of the the probability 



density function ( PDF ) description of WT theory, we explain the coexistence and mutual 



interactions between solitons and the weakly nonlinear random wave background in the 



form of a wave turbulence life cycle ( WTLC I 



Keywords: Nonlinear optics, liquid crystals, turbulence, integrability, solitons 



Contents 

1 Introduction [4] 

2 The Experimental Setup of the One-Dimensional Optical System [8] 



* Corresponding author 

Email addresses: jason.laurie@ens-lyon.fr (Jason Laurie), umberto.bortolozzoOunice.fr 
(Umberto Bortolozzo), s.v.nazarenko@warwick.ac.uk (Sergey Nazarcnko), 
stefania.residori@inln.cnrs.fr (Stefania Residori) 

Preprint submitted to Elsevier October 24, 2011 



3 The Relation to Previous Studies of Optical Solitons [9] 



The Theoretical Model of the |OWT| Experiment ]W\ 



4.1 The Long- Wave Regime [TT] 

4.2 The Short- Wave Regime [T2J 

4.3 The Nonlinearity Parameter Q2] 

4.4 The Hamiltonian Formulation Q2] 

5 Wave Turbulence Theory 1151 

5.1 The Statistical Setup [15] 

5.2 The Canonical Transformation [TH] 

5.3 The Weak Nonlinearity Expansion [lj|] 

5.4 Equation for the Generating Functional [21] 

5.5 Moments of the Wave Intensity [22] 

5.6 The One-Mode Amplitude Probability Density Function [22] 

5.7 The Kinetic Equation [23] 

5.8 The Fj0rtoft Argument [22] 

5.9 The Zakharov Transform and the Kinetic Equation Solutions [25] 

5.10 Locality of the Kolmogorov- Zakharov Solutions [29] 

5.11 Logarithmic Correction to the Direct Energy Spectrum [30] 

5.12 Linear and Nonlinear Transfer Times and The Critical Balance Regime . [3JJ 

5.13 The Differential Approximation for the Kinetic Equation [32] 

5.14 Modulational Instability and Solitons in the Long- Wave Equation [34] 

6 The Numerical Method 

7 Results 

7.1 The Long- Wave Equation [UJ 

7.1.1 The Decaying Inverse Cascade with Condensation Hi] 

7.1.2 Forced and Dissipated Simulations [5S] 

7.2 The Short- Wave Equation 1JJ 



8 Conclusions 
Appendices 

A The Canonical Transformation 

B Phase and Amplitude Averaging 1721 

C The Zakharov Transform [75] 

D Expansion of the Long- Wave Six- Wave Interaction Co- 
efficient [76] 

E Locality of the Kinetic Equation Collision Integral 1771 

F Derivation of the Differential Approximation Model 1791 



2 



G The Bogoliubov Dispersion Relation 1801 

H Non-Dimensionalization 1811 

I The Intensity Spectrum [82] 

Bibliography [83] 



3 



1. Introduction 



One-dimensional ( ID ) optical wave turbulence ( OWT I is an extremely interesting 



physical phenomenon whose importance arises from its intrinsic overlap with several 
strategic research areas. This interdisciplinary nature allows for the application of non- 
conventional approaches to familiar facts and routes. These areas include wave turbulence 



( WT ) , Bose-Einstein condensate ( BEC I and lasing, integrable systems and solitons, and, 



on a more fundamental level, general turbulence, nonlinear optics, equilibrium and non- 
equilibrium statistical mechanics. A hierarchical diagram showing these areas and their 
links to |lD|[OWT| is shown in figure [T] 
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Figure 1: Interconnections of |1D|QWT| with other research areas. 



The main aspects and phenomena in |1D||QWT| which will be at the focus of the 
present review include: 

• The inverse transfer of wave action (from short to long wave lengths), its relation 
to the BEC of light, and to WT| as an active medium for lasing; 



• The role of turbulent cascades (fluxes) versus thermodynamic potentials (temper- 
ature and chemical potential) in 1D||QWT and in WT in general; 



• The proximity to integrability, modulational instability ( MI ) and the generation of 
solitons; 

• The coexistence, interactions and mutual transformations of random waves and 



solitons: this constitutes the wave turbulence life cycle (WTLCl and the evolution 
towards a final single soliton state. 



Since |WT| conc epts are central for our review, we will begin by giving a brief in- 
troduction to WT WT can be generally defined as a random set of interacting waves 
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with a wide range of wave lengths. |WT| theory has been applied to several physical 
systems including water surface gravity and capillary waves in oceans [H [^1 El SI El IH1 
El El El HOI El H21 H31 [HI HH US]; internal, inertial, and Rossby waves in atmospheres 
and oceans [13 HH1 CHI HOI HH HU 113 HH US] ; Alfven waves in solar wind and interstellar 
turbulence [26l ETJ Ell EH EOl EU [32l EH S3 EH EH (37]; Kelvin wa ves on quantized 
vortex lines in superfluid helium ESI E3 [4Ql SH [42l [43l [44] ; waves i n|BEC[ > and non- 
linear optics [3S1 SHI H71 SSI US] ; waves in fusion plasmas EkH ED EH E|; and waves on 
vibrating, elastic plates |53) . A thorough and detailed list of examples where the WT 



approach has been applied, from quantum to astrophysical scales, can be found in recent 
book [54] . The most developed part of the WT theory assumes that waves have random 
phases and amplitudes and that their interactions are weakly nonlinear El .54], in which 
case a natural asymptotic closure arises for the statistical description of |WT| The most 
familiar outcomes of such a closure is the kinetic equation ( KE ) for the wave action spec- 
trum and its stationary solutions describing energy and wave action cascades through 



scales called Kolmogorov-Zakharov ( KZ I spectra El E3] • It is the similarity of the KZ 



spectra to the Kolmogorov energy cascade spectrum in classical three-dimensional ( 3D ) 
Navier-Stokes turbulence that allows one to classify |WT| as turbulence. On the other 
hand, in most physical applications, besides weakly nonlinear random waves, there are 
also strongly nonlinear coherent structures which, even when not energetically dominant, 
interact with the random wave component, i.e. by random-to-coherent and coherent-to- 
random transformations, which may provide a route to the turbulence sink via wave 
breaking or wave collapses. In other words, together with the random weakly nonlinear 
waves, the strongly nonlinear coherent structures are also a fundamentally important 
part of the|WTLC]E3]. 



Let us now discuss realizations of WT in optical systems. Very briefly, we define 



OWT as the WT of light. As such, OWT is a subject within the more general nonlinear 



optics field, dealing with situations involving the propagation of light in nonlinear media 
which is fully or partially random. |OWT| is a niche area within the nonlinear optics field, 
and excludes are large section of the field, including systems fully dominated by strongly 
coherent structures, such as solitons. The term turbulence, is even more relevant for 



OWT because of the similarities between the nonlinear light behavior to fluid dynamics, 

and shock waves [57 . Although there have been 



such as vortex-like solutions 1551 



numerous theoretical and numerical studies of OWT [48, 58] EH E01 EI] > there have been 
few experimental observations to date [35] ■ OWT was theoretically predicted to exhibit 



dual cascade properties when two conserved quantities cascade to opposite regions of 
wave number space |48j . This is analogous to two-dimensional ( 2D I turbulence, where 
we observe an inverse cascade of energy and a direct cascade of enstrophy [611 162] . In 
the context of [OWTj energy cascades to high wave numbers, while wave action cascades 
towards low wave numbers [48 , 58 , 59 , 60J . An interesting property of OWT is the inverse 
cascade of wave action which in the optical context implies the condensation of photons 



the optical analogue of BEC 



It is the BEC processes that make OWT an attractive and important setup to study. 
Experimental implementation of [BEC in alkali atoms was first achieved in 1995, and 
subsequently awarded the 2001 Nobel prize E21E3]- This work involved developing a 
sophisticated cooling technique to micro-Kelvin temperatures, in order to make the de 
Broglie wave length exceed the average inter-particle distance. This is known as the [BEC] 
condition. Photons were actually the first bosons introduced by Bose in 1924 [65 , and 
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the |BEC] condition is easily satisfied by light at room temperature. However, there was 
the belief that optical [BEC| would be impossible, because of the fundamental difference 
between atoms, whose numbers are conserved, and photons which can be randomly 
emitted and absorbed. However, there exist situations where light is neither emitted 
nor absorbed, e.g. light in an optical cavity, reflected back and forth by mirrors [66) . or 
by light freely propagating through a transparent medium. In the latter, the movement 
of photons to different energy states (specifically the lowest one corresponding to BEC ) 
can be achieved by nonlinear wave interactions. The mechanism for these nonlinear 
wave interactions is provided by the Kerr effect which permits wave mixing. Moreover, 
the nonlinear interactions are crucial for the BEC of light, because no condensation is 



possible in non-interacting [ID] and [2D] Bose systems. 

When the nonlinearity of the system is weak, |OWT| can be described by weak |WT| 
theory [3] [53] , with the prediction of two KZ states in a dual cascade system. One aspect 
of |OWT| is that the nonlinearity of the system is predicted to grow in the inverse cascade 
with the progression of wave action towards large scales. This will eventually lead to a 
violation of the weak nonlinearity assumption of |WT| theory. The high nonlinearity at low 
wave numbers will lead to the formation of coherent structures [48 , 59 , 45 , 67, 68, 69J. In 
|OWT| this corresponds to the formation of solitons and collapses for focusing nonlinearity 
[IH], or to a quasi- uniform condensate and vortices in the de- focusing case [45] . 

Experimentally, OWT| is produced by propagating light through a nonlinear medium, 
e.g. a photo-refractive crystal. However, the nonlinearity is typically very weak and it is 
a challenge to make it overpower the dissipation. This is the main obstacle regarding the 
photon condensation setup in a |2D| optical Fabry-Perot cavity, theoretically suggested in 
[55] . Recently the problem of dissipation in the Fabry-Perot cavity has been elegantly 
approached by using a very short cavity filled with a drop of an active medium, thus 
providing gain and allowing photon thermalization with subsequent observation of con- 
densation [70) . However, this effect encompasses a gain mechanism and, hence, is outside 
the |OWT| regime. 

This brings us to discussing an exceptional role played by |1D| optical systems - |1D| 



OWT Firstly, it is in ID that the first ever OWT experiment was implemented [41?] . 
The key feature in our setup is to trade one spatial dimension for a time axis. Namely, 
we consider a time-independent 2D light field where the principal direction of the light 
propagation acts as an effective time. This allows us to use a nematic liquid crystal (LC I, 
which provides a high level of tunable optical nonlinearity |711 172] . The slow relaxation 
time of the re-orientational dynamics of the |LC| molecules is not a restriction of our 
setup because the system is steady in time. Similar experiments were first reported in 
[73], where a beam propagating inside a nematic layer undergoes a strong self- focusing 
effect followed by filamentation, soliton formation and an increase in light intensity. 
Recently, a renewed interest in the same setup has led to further studies on optical 
regime [71 [75] [75]. 



solitons and the MI 



However, the previous experiments used a high 
input intensity, implying a strong nonlinearity of the system, and therefore the soliton 
condensate appears immediately, bypassing the WT regime. In our experiment, we set 



up an initial condition of weakly nonlinear waves situated at high wave numbers from a 
laser beam. We randomize the phase of the beam, so that we produce a wave field as 
close to a random phase and amplitude (RPA) wave field as possible. The nonlinearity 



of the system is provided by the |LC| controlled by a voltage applied across t he|LC| cell. 
This provides the means for nonlinear wave mixing via the Kerr effect. The |LC| we use 
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is of a focusing type, causing any condensate that forms to become unstable and the 
formation of solitons to occur via IMI1 

Secondly, from the theoretical point of view |1D||QWT| is very interesting because it 
represents a system close to an integrable one, namely the ID nonlinear Schrodinger 
Thus it inherits many features of the integrable model, e.g. the 

On the other hand, de- 



equation (NLSE). 



significant role of solitons undergoing nearly elastic collisions 
viations from the integrability are important, because they upset the time recursions of 
the integrable system thereby leading to turbulent cascades of energy and wave action 
through scales. The process responsible for such cascades is a six-wave resonant inter- 
action (wave mixing). Another example of a nearly integrable ID six- wave system can 
be found in superfluid turbulence - it is the |WT| of Kelvin waves on quantized vortex 
lines [311 [77J EH H31 HI] (even though non-local interactions make the six-wave process 
effectively a four- wave one in this case). Some properties of the six- wave systems are 
shared with four- wave ID particularly WT in the Majda-McLaughlin-Tabak (MMT) 



model reviewed in Physical Reports by Zakharov et al. 69 . For example, both four- 
wave and the six-wave systems are dual cascade systems, and in both systems solitons 
(or quasi-solitons) play a significant role in WTLC There are also significant differences 
between these two types of systems. Notably, pure |KZ| solutions appear to be much 
less important for the six- wave optical systems than for the four- wave |MMT| model - in- 
stead the spectra have a thermal component which is dominant over the flux component. 
Moreover, the number of solitons in ID OWT decrease in time due to soliton mergers, 
so that asymptotically there is only a single strong soliton left in the system. 

Similar behavior was extensively theoretically studied in various settings for non- 
integrable Hamiltonian systems starting with the paper by Zakharov et al. [75], and 
then subsequently in [7H1 [HOI Ell IMl IH31 IB11 HS1 [13 HB] ■ The final state, with a single 
soliton and small scale noise, was interpreted as a statistical attractor, and an analogy was 
pointed out to the over-saturated vapor system, where the solitons are similar to droplets 
and the random waves behave as molecules [88] . Indeed, small droplets evaporate whilst 
large droplets gain in size from free molecules, resulting in a decrease in the number of 
droplets. On the other hand, in the |1D|[0WT| context, the remaining strong soliton is 
actually a narrow coherent beam of light. This allows us to interpret the WT evolution 
leading to the formation of such a beam as a lasing process. Here, the role of an active 
medium where the initial energy is stored is played by the weakly nonlinear random wave 
component, and the major mechanism for channeling this energy to the coherent beam 
is provided by the |WT| inverse cascade. For this reason, we can call such a system a 
WT\ laser. It is quite possible that the described WT lasing mechanism is responsible 



for spontaneous formation of coherent beams in stars or molecular clouds, although it 
would be premature to make any definite claims about this at this point. 

Generally, in spite of recent advances, the study of |1D|[Q~WT| is far from being com- 
plete. The present review provides a report on the current state of this area describing 
not only what we have managed to learn and explain so far, but also the results which 
we do not know yet how to explain, discussing the existing theory and the gaps within 
it that are yet to be filled. We compare the experimental observations with the pre- 
dictions of |WT| theory and independently juxtapose our findings with direct numerical 
simulation of the governing equations. In particular, we describe some puzzles related to 
the wave action spectra obtained in the experiment and in the numerical experiments. 
Furthermore, we will describe the recent extensions of the weak WT| theory onto the 
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wave probability density function ( PDF I , which marks the beginning of developing a 
formalism for describing |WT| intcrmittency and the role of the coherent structures. On 
the other hand, a theory for the |WTLC| incorporating interacting random waves and 
coherent structures/solitons remains largely incomplete. 

2. The Experimental Setup of the One-Dimensional Optical System 

The |1D| optical system has been designed to meet the major requirements of |OWT| 
Especially important were the careful calibrations that have been taken to fulfill the 
balance between low dissipation and low nonlinearity. Indeed, the main experimental 
challenge to observe the |WT| regime is that of keeping the nonlinearity weak enough to 
let the |WT| regime develop and, at the same time, high enough to make it overpower 
the dissipation. Our setup is based on a nematic liquid crystal layer in which a laminar 
shaped beam propagates. |LCfc are particularly suitable for the observation of the |WT| 
regime because of their well known optical properties, such as their high and tunable 
nonlinearity, transparency (slow absorption) over a wide range of optical wave lengths, 
the realization of large cells and the possibility to drive them with low voltage externally 
applied fields [75]. 

The experimental apparatus is shown in figure [2^i. It consists of a |LC| cell, inside 
which a laminar shaped beam propagates. The input beam is prepared in such a way, 



as to have an initial condition of weakly nonlinear random waves. The LC cell is made 
by sandwiching a nematic layer, (E48), of thickness d = 50 /im, between two 20 x 30 
mm 2 , glass windows and is schematically depicted in figure ^p. On the interior, the 
glass walls are coated with indium-tin-oxide transparent electrodes. We have pre-treated 
the indium-tin-oxide surfaces with polyvinyl-alcohol, polymerized and then rubbed, in 
order to align all the molecules parallel to the confining walls. When a voltage is applied 
across the cell, |LC| molecules tend to orientate in such a way as to become parallel to 
the direction of the electric field. By applying a 1 kHz electric field with a rms voltage 
of Vq = 2.5 V we preset the molecular director to an average tilt angle O. 

The |LC layer behaves as a positive uni-axial medium, where nn — n z = 1.7 is the 



extraordinary refractive index and nj_ = 1.5, the ordinary refractive index |72| . The LC 
molecules tend to align along the applied field and the refractive index, n(8), follows 
the distribution of the tilt angle 6. When a linearly polarized beam is injected into the 
cell, the |LC| molecules orientate towards the direction of the incoming beam polarization. 
The input light originates from a diode pumped, solid state laser, with Ao = 473 nro, 
polarized along y and shaped as a thin laminar Gaussian beam of 30 \im thickness. 
The beam evolution inside the cell is monitored with an optical microscope and a CCD 
camera. The light intensity is kept very low, with an input intensity of I — 30 fiW/ cm 
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to ensure the weakly nonlinear regime. A spatial light modulator (SLM), at the entrance 
of the cell is used to produce suitable intensity masks for injecting random phased fields 
with large wave numbers. 

Figure [3] depicts how the input light beam is prepared before entering the |LC| cell. 
The beam is expanded and collimated through the spatial light filter shown in figure [3] 
The objective, OB, focuses the light into the 20 /im pinhole, PH, the lens Li collimates 
the beam with a waist of 18 mm. After that, the light passes through the |SLM} which is 
a LCD screen working in transmission with a resolution of 800 x 600, with 8 bits pixels, of 
size 14 nm. Each pixel is controlled through a personal computer PC, to ensure that the 




Figure 2: a) Sketch of the experimental setup: a laminar shaped input beam propagates inside the 
|LC| layer; random phase modulations are imposed at the entrance of the cell by means of a SLM. b) 
Schematic of the ILCI cell. 



outgoing light is intensity modulated. In our case we use a cosinus modulation having a 
colored noise envelope. The lenses, L3 and L4, are used to focus the image from the LCD 



screen at the entrance of the LC cell. The half wave-plate, W, together with the polarizer, 
P, are used to control the intensity and the polarization, which is linear along the y-axis. 
The circular aperture is inserted in the focal plane to filter out the diffraction given by 
the pixelization of the |SLM| and the diffuser, PH, is inserted to spatially randomize the 
phase of the light. In order to inject the light inside the LC cell, we use a cylindrical 



lens, L4, close to the entrance of the LC layer 
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Figure 3: Detailed representation of the experimental setup, showing the initialization of the input laser 
beam. 



3. The Relation to Previous Studies of Optical Solitons 



Optical reorientation of the [LC] occurs under the action of the light itself, with the |LC 
molecules tending to align along the direction of the laser beam polarization, and giving 
rise to a Kerr effect that produces a self- focusing nonlinearity [7T] . This effect was shown 
to lead to laser beam filamentation j73j and, recently, has been exploited to demonstrate 
the stable propagation of spatial optical solitons, also called nematicons, inside nematic 



LC cells [71]. For the same type of system, MI has also been reported [751 [75]. However, 



the previous experiments used a high input intensity, implying a strong nonlinearity of 
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the system, and therefore filamentation and solitons appear immediately, bypassing the 
|WT| regime. 

Solitons are well-known and widely studied in nonlinear optics |89j . where they are 
understood as light pulses that maintain their shape unaltered, while propagating in a 
nonlinear medium, and where nonlinearity means a change in the refractive index that 
is induced by the intensity of the light itself. In this context solitons are classified as 
being cither spatial or temporal, whether the self-confinement of the light beam occurs in 
space or in time during their propagation. The nonlinearity of the medium corresponds, 
respectively, to a self-focusing or a self phase modulation effect. 

While |WT| implies the presence of many random waves interacting with low nonlinear- 
ity, previous work conducted in the nonlinear optics field was mainly aimed at realizing 
self-confined beams in a different regime, where the nonlinearity was relatively high and 
the medium seeded by a single pulse of Gaussian shape. A spatial soliton is therefore 
obtained by imposing a tightly focused beam as an initial condition and, then by letting 
the beam propagate inside a medium where the self-focusing nonlinearity compensates 
the transverse beam widening by diffraction. These types of spatial solitons have been 
observed in a number of diverse optical media, such as, photo-refractive crystals, atomic 
vapors, and semiconductor wave guides [301 HI]- O n the other hand, temporal solitons 
correspond to cases when the nonlinearity compensates the temporal broadening of a 
light pulse due to the natural dispersion of the traversed medium 72], and have been 
observed in optical fibers. 

In both cases, the theoretical approach to spatial and the temporal solitons are based 
on the |NLSE| We therefore expect the generic behavior of the system to have an |OWT| 
regime existing before the formation of solitons, then with an inverse cascade, the de- 
velopment of soliton turbulence and subsequently, a final single soliton as a statistical 
attractor of the system |78j . In the spatial case, the light propagation direction, usually 
denoted as z, plays the role of time, therefore the wave dispersion relation is of the type 
to = k 2 , whereas in the temporal case the time derivative is in the second order disper- 
sion term, hence the dispersion relation is of the type k = uj 2 . Theoretical predictions 
of |QWT| regimes should, therefore, be different in both cases. Recently, experiments 



achieving WT regimes in optical fibers have been devised |92j and the wave thermaliza- 
tion phenomena has been reported for the same type of system |93j . However, up to 
now only in the |LC| experiment has there been a genuine |WT| regime, with an inverse 
cascade and the spontaneous emergence of spatial solitons out of random waves, been 
demonstrated [4"9"] . 

4. The Theoretical Model of the |OWT| Experiment 

Theoretically, the experimental setup can be modeled by an evolution equation for 
the input beam, coupled to a relaxation equation for the |LC| dynamics given by 



2iq jz + di + (la) 

d 2 a 1 e n 2 2 

^-q a+ ^ =°< (lb) 

where ip(x, z) is the complex amplitude of the input beam propagating along the time 



axis z; x is the coordinate across the beam; a(x, z) is the LC reorientation angle; n a 
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n e — n is the birefringence of the |LC| fc is the optical wave number; e is the vac uum 
permittivity; and — J nK/2Ae(d/Vo) is the electrical coherence length of the 



LC 



with K being the elastic constant, q 2 = k 2 (n 2 + n 2 a /$) and Ae is the dielectric anisotropy. 
Note that It fixes the typical dissipation scale, limiting the extent of the inertial range 
in which the OWT cascade develops. In other contexts, such a spatial diffusion of the 
molecular deformation has been denoted as a non-local effect, see [74l [75l [76] . In our 
experiment, when Vq = 2.5 V, we have that 1^ — 9 /im. By considering that a typical 
value of K is of the order ~ 10 pN, we can derive a typical dissipation length scale of 
the order ~ 10 /im. 

The evolution equations, ([!]) for the complex amplitude of the input beam, ip(x,z), 
to the LC reorientation angle, a(x, z) can be re-written as a single equation for field 
ip(x,z). This can be achieved by formally inverting the operator on a(x, z) in equation 



( lb ) . Equation ( lb I implies that 



a(x, z) 




(2) 



Substituting this expression into equation (la), we eliminate the dependence on variable 
a(x,z). This gives 




(3) 



Equation ^ is a single equation modeling the evolution of the complex amplitude, 
ip(x,z). We can further simplify equation ([3| by considering the system in two limits 
of wave number k: kl^ <C 1 and 1 kl^, that we call the long- and short-wave limits 
respectively. These limits enable the expansion of the nonlinear operator in power of 
l^d/dx. Our experimental system is well described by the long-wave limit. The limita- 
tions imposed by the dissipation of the |LC| in the current experimental setup prevents 
the implementation in the short-wave regime. However, for the completeness of our de- 
scription and for the possibilities in the modification of the experimental setup in the 
future for the short-wave regime, we will continue to investigate this limit theoretically 
and numerically. 

4-1- The Long-Wave Regime 

The long-wave approximation to equation (|3| corresponds to the wave length of the 
light field in the LC cell, A oc 1/fc, being greater than the electrical coherence length of 
the LC l^. In physical space, this limit corresponds to l^d/dx -C 1, which permits the 
expansion of the nonlinear operator of equation ^ as 



(4) 




Taking the leading order of this expansion yields 



(5) 



where, for clarity, we have introduced a reference light intensity: 



2K 
e nil\kl ' 



(6) 



Equation ([5| is the ID focusing NLSE As is well-known, the 1D||NLSE is an integrable 
system, solvable with the aid of the inverse scattering transf orm [95 , and characterized 
by soliton s [[| Un fortunately, this would be a poor model for OWT as the integrability 
of the 1D||NLSE implies that wave turbulent interactions are not possible. To overcome 
this, we must consider the sub-leading contribution in expansion Q. This extra nonlinear 
term acts as a correction breaking the integrability of the system. The resulting equation 



is given as 



2iq 



dip 



d 2 ip 
dx 2 



2lU 2 



V#l 2 



(7) 



We refer to equation Q as the long- wave equation ( LWE I . As a consequence of expansion 
Q, the additional nonlinear term is smaller than the leading nonlinear term. Moreover, 
the additional, sub-leading, nonlinear contribution is smaller than the linear term in 
the weak nonlinearity setting. Therefore, although integrability is lost, the system will 
remain close to the integrable one described by |5j). As a result, we expect soliton-like 
solutions close to the exact solutions of the ID |NLSE| (j5J). On the other hand, exact 
soliton solutions of equation ([5]), do not change shape, and have the ability to pass 
through one another unchanged. We can expect that in the |LWE equation Q, we will 
observe similar soliton solutions, but the non-integrability will allow solitons to interact 
with one another, and with the weakly nonlinear random wave background. 

4-. 2. The Short- Wave Regime 

In the opposite limit of equation (|3j), when 1 -C l 2 d 2 /dx 2 1 the nonlinear operator of 
equation Q can be represented in terms of a Taylor expansion of negative powers of the 
spatial derivative: 




dx 2 



I 2 



dx 2 



(8) 



It is sufficient for us to approximate the nonlinear operator of equation ([3| with just 
the leading order term in expansion ([8]), as integrability of the equation is not an issue. 
Therefore, we get an equation of the form: 



2iq 



dtp 
dz 



d 2 ip 
dx 2 



2lU 2 



dx 2 



(9) 



We call equation ^ the short-wave equation (SWE) 



two dynamical equations for the complex wave amplitude tp(x, z) for ID OWT 



Ultimately, we have presented 
in two 



limits of wave number space. Both of these systems can be expressed in a Hamiltonian 
formulation, that will be utilized by |WT| theory in the weakly nonlinear regime. 



The term soliton is sometimes reserved for solitary waves with special properties arising from inte- 
grability, such as the ability to pass through one another without change in shape or velocity, as is the 
case for the |lD|NLSE| Hereafter, we will use the term soliton more broadly, including solitary waves in 
non-integrable systems which can change their states upon mutual collisions. 
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4-3. The Nonlinearity Parameter 

To verify that experimental and numerical simulations are in a weakly nonlinear 
regime, we can quantify the linearity and nonlinearity within the system with the intro- 
duction of a nonlinear parameter, J, which is determined by the ratio of the linear term 
to the nonlinear term within the dynamical equations. 

The nonlinear parameter from the LWE| Q is defined as 



, _ 2^ 2 k 2 l 2 



(10) 



This is derived from the ratio of the linear term and the first of the two nonlin ear ter ms. 
Here, / = (\tp(x, z)| 2 ) is the average value of the light intensity. Similarly, the SWE ^ 
yields a nonlinearity parameter of 



J s = 



24> 2 k 4 lj 



(11) 



Calculation of JL and J s act as a verification of the weak nonlinear assumption of 



WT especially in the context of the experimental setup where unknown quantities are 



more difficult to measure. 



4-4- The Hamiltonian Formulation 

Both equations |7]) and Q can be written in terms of a Hamiltonian system of the 

dtp SH 



form 



2iq 



dz Stp* 



(12) 



For the LWE[ the Hamiltonian is given as 



njL nj i njL 

tt — til T rt.4 , 





dtp 


2 


1 




dx 




4^ 2 



IV>| 4 fdW 



;2 



V dx 



dx. 



(13a) 



In the nonlinear energy term %4, the term quartic with respect to tf>, we have added 
a superscript L to denote that this quartic term corresponds to the LWE| Q. This is 
because the Hamiltonian of the |LWE| and |S WE| only differ in the expression H4. For the 
|SWE[ the Hamiltonian is given by 



n 2 



dtp 

dx 



4 ! 

2 



1 1 

4Zp2 ^ dx- 1 



dx. 



(13b) 



In both the LWE and the SWE, the linear, (quadratic), energy H.2 is identical. The 
Hamiltonians ( 13 ) coincide with the total energy of the systems and are conserved by 
their respective dynamics (H — const). Moreover, both the LWE and SWE contain an 
additional invariant, the wave action J\f defined as 



Af = J \tp\ 2 dx. 
13 



(14) 



Conservation of A/" is a consequence of the U(l) gauge symmetry or invariance with 
respect to a phase shift: ip(x, z) — > tf)(x, z) exp (i<f>). 

By expressing the Hamiltonian in terms of its Fourier representation 



(15) 



here k G 
as 



the general Hamiltonian structure for Hamiltonians ( 13 ) can be represented 
^w fe a k a£ + - T 3,4 axa 2 a%al 8\% (16) 

k 1,2,3,4 



H 



-k 4 ) 



where u>k = k 2 is the linear frequenc}|^|of non-interacting waves, \ = <5(ki+k 2 — k 3 - 
is a Kronecker delta function, and T 3 ' 4 = T(ki, k2, k3, k 4 ) is the nonlinear interaction 
coefficient. Note that we use bold symbol k for the wave number to emphasize that it 
can be either positive or negative, while k is reserved specifically for the wave vector 
length, k = |k|. 

By symmetry arguments, the interaction coefficient should not change under the 
permutations ki -h- k2, k 3 4-> k 4 . Furthermore, the Hamiltonian (16) represents the 
total energy of the system and is therefore a real quantity. This property implies extra 
symmetries of the interaction coefficient: 



rpl.2 _ rp2,l _ rpl.2 _ fm3,4\< 

1 3A — J 3.4 — J 4,3 — V J 1,2,1 



(17) 



For the LWE Hamiltonian (|13a|), the interaction coefficient is defined as follows, 



1.2 



3.4 



J- 1 A T J- 1 



-,1,2 

L 3,4 

1 



2^1,2 

3,4 



-i- (kik4 

2^ 2 



k 2 k 3 + k 2 k 4 + kik 3 - 2k 3 k 4 - 2kik 2 



(18a) 



ID 



We have denoted the two contributions to L T 3 1 4 , from both nonlinear terms in (|7j), as 

1 T^'2 and 2 T\^ , where the first arises from the usual cubic nonlinearity seen in the 
focusing |NLSE[ and the second from the sub- leading correction. 
Similarly, the [SWE] yields the following interaction coefficient, 



3,4 



kik 3 k 9 k 



2*3 



kik 4 k 5 k, 



2*4 



kik 2 



k 3 k 4 



(19a) 



In terms of the wave amplitude variables a(k) 
Fourier analogue of equation p2| : 

da(h, z) 

dz 



the Hamiltonian system ( 16 1 satisfies the 



2iq 



sn 



5a*(\t,z) 



(20) 



It is with equation (20) that the formulation of WT theory is applied. In the next 
Section, we will give a brief mathematical description of |WT| theory, and outline the 
assumptions on the wave field that is required to apply such an approach. 



2 Indeed, is the frequency with respect to the time variable which is related to the distance z as 
t = z/2q. 
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5. Wave Turbulence Theory 



Historically, there have been several approaches taken in the formulation of |WT| 
theory, initially this was achieved using a diagrammatic technique [M] EE], then 
by cumulants expansions [991 H00| H01| 11021 1103j , and finally by the use of the random 
phase approximation 3, 2 . However, the random phase approximation is insufficient for a 
rigorous treatment of |WT| theory - the wave amplitudes must also be random independent 
variables. A formulation of |WT| theory, where the phases and the amplitudes are random, 
and yet the wave fields do not necessarily have to be Gaussian, was recently developed 
in (1041 11051 11061 11071 [54] . a strategy which we will be followed in this review. This 



approach allows for the derivation of not only the |KE| but also an equation for the [PDF] 
and for the ability to consider non-Gaussian wave fields. These advantages provide a 
formulation for the study of the underlying |WT| assumptions, as well as for description 
of the WT intermittency and to extending the theory to incorporate the presence of 
coherent structures. 



5.1. The Statistical Setup 

The general multi-dimensional formulation for WT theory can be found in |54) , but for 
the treatment of |lD|OWT| and for reasons of clarity, it will be sufficient to limit ourselves 



to presenting the WT description in ID Let us consider a [ID] wave field, a(x,z), in a 
domain which is periodic in the ^-direction with period L, and let the Fourier transform 
of this field be represented by the Fourier amplitudes <z;(t) = a(k/, z), with wave number, 
k; = 2ttI/L, I 6 1. Recall that the propagating distance z in our system plays a role of 
"time", and consider the amplitude-phase decomposition a;(t) = A l (z)tfji(z), such that 
Ai is a real positive amplitude and ipi is a phase factor that takes values on the unit 
circle in the complex plane. 



We will begi n by fo llowing [TUSl QUO DM HMj and define a |RPA| field. We say a wave 
field a(x,z) is a RPA field, if all the amplitudes Ai{z) and the phase factors ipi(z) are 



independent random variables and all tp^s are uniformly distributed on the unit circle 
on the complex plane. We remark, that the |RPA| property does not require us to fix 
the shape of the amplitude [PDF] and therefore, we can deal with strongly non-Gaussian 
wave fields. This will be useful for the description of WT intermittency. By Constructing 



a generating functional (GF), and subsequently deriving an equation for its evolution, 



we can systematically derive dynamical equations for all the statistical quantities of the 
wave system. 

Here, we will only concentrate on the one-mode statistics of the system. Therefore, 



we will only need to consider the one-mode amplitude GF defined as 



2k (Ai, 



(21) 



where Jk — \a^\ 2 for wave number k. Besides from deriving all the statistical moments of 
the wave field from the |GF[ we can also obtain the one-mode |PDF| for wave intensities 
Jk by taking the inverse Laplace transform of the one-mode amplitude |GF[ 



^k(sk) 



1 

2ni 



/ e -(-) A >^ Zk(Ak) d\ v 

J — zoo+c 



(22) 



where c is a constant greater than the real part of all singularities. 
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So far we have presented the necessary statistical properties of the underlying wave 
field for the application of |WT| theory. We now present the application of the theory to 
the |1D||QWT| Hamiltonian system presented in Section [4] 

5.2. The Canonical Transformation 

Nonlinear wave interactions can be classified by the lowest order of resonance interac- 
tions they undergo. For anJVf>M wave scattering process, these resonance conditions 
are defined as 

kH h k w = k w+ i H h kjy+M, (23a) 

£*>l H h u N = u>n+i H h u N+M , (23b) 

where k^ is the wave number and u>i — w(kj) is the frequency of wave i. 

Of course, apart for satisfying the resonant conditions, the respective type of the non- 
linear coupling must be present. For example, the |2D| and |3D||NLSE| have the dispersion 
relation ui^ = k 2 which can satisfy the three-wave lf)2 resonance conditions, but the 
three-wave nonlinear coupling is zero. On the other hand, for the gravity water waves 
there is a 1 H 2 wave interaction Hamiltonian (when written in terms of the natural 
variables - height and velocity potential), but the wave linear frequency Wfc = y/gk does 
not allow for 1 o 2 resonances [50]. As a consequence, the lowest order resonant pro- 



cesses in all of these cases are four- wave (2f> 2). For ID OWT like in the NLSE the 
frequency of the linear propagating waves is given by 

w(k) = k 2 , (24) 



In ID dispersion relations of the form w(k) oc k a with a > 1 cannot satisfy the four-wave 
resonance condition 

k 1 +k 2 = k 3 + k 4 , (25a) 

w(ki)+w(k 2 ) =u;(k 3 )+u;(k4). (25b) 

( |25[ ). This can be understood by a simple graphical proof presented in figure |i]|5ni E3I- In 
figure [4] we observe the red dashed curve representing the dispersion relation to^ = Ck a 
with a > 1. At two locations along this curve (at k = ki and k = k 3 ), two further 
dispersion curves (the green and blue solid lines) are produced: with their minima at 
points (k, w) = (ki,wi) and at ^3,^3) respectively. These subsequent lines represent 
the wave frequencies of u>i + u; 2 and W3 + W4, where ki and k3 are now fixed, with k 2 
varying along the green solid line and k4 varying along the blue line. If the green and blue 



lines intersect, it will be when the four-wave resonance condition (25 1 is satisfied and will 
occur at the point (k, w) = (ki +k2, ui\ +0J2) = 0^-3 +k4, ui 3 +o; 4 ). In figure|3j we observe 
that such an intersection occurs only once, and it can be clearly seen that ki = k4 and 
k 2 = k3 must hold. This corresponds to a trivial pairing of wave numbers, which will not 
provide any nonlinear energy exchange between modes. As a consequence, resonant four- 
wave interactions are absent in the system. There are no five-wave interactions either 
because the U(l) symmetry prohibits the presence of odd orders in the wave amplitude in 
the interaction Hamiltonian. In situations such as these, there exists a weakly nonlinear 



canonical transformation ( CT ) that allows us to change to new canonical variables such 



that the leading interaction Hamiltonian is of order six. 



1G 



Figure 4: We plot a graphical representation of the four- wave resonance condition |54| . The four- wave 
resonance condition is satisfied at points where the green and blue lines intersect, (shown by the black 
dot). However, for dispersion relations ui^ <x k Q , with a > 1, there can only be one intersection, 
corresponding to the trivial wave resonance: ki = k4 and k2 = k3. 



Such a |CT| is similar to the algorithm in the construction of the corrected wave 



action for the perturbed integrable systems in KolmogorovArnoldMoser ( KAM ) theory. 
The latter represents a recursive procedure eliminating the lower-order interaction terms 
one by one at each of the recursive steps, which is only possible when there are no 
resonances at that respective order. In |WT| theory, such a recursion is "incomplete" - it 
contains a finite number of recursive steps until the later steps are prevented by the lowest 
order wave resonances. The |CT| procedure for eliminating the non-resonant interactions 
in WT theory is explained in [3J, and the most prominent example was given for the 
system of the gravity water waves, where it was used to eliminate the non-resonant cubic 
Hamiltonian (see also |109j where some minor mistakes made for surface gravity waves 
were corrected). 

A similar strategy was recently applied to eliminate non-resonant fourth-order inter- 
actions in the context of Kelvin waves in superfluid turbulence [1101 144) and in nonlinear 



optics [49] . The details of the CT for our optical system can be found in Appendix 
The result of the CT is the representation of our system in a new canonical variable 
Ck with the elimination of the quartic contribution H4. This however results in the 
appearance of a sextic contribution 'Hg : 

H = ^W fc C k C k + — ^ W 4,5,6 ^4,5,6 CiC 2 C 3 cl44, (26) 
k 1,2,3,4,5,6 
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Figure 5: An illustration to show the non-resonant four- wave interaction, T 1 '^ , and the resonant six- wave 
interaction, W^'g'g, after the CT The six- wave (sextet) interaction term is a contribution arising from 
two coupled four- wave (quartet) interactions via a virtual wave (dashed line). 



where the explicit formula for W^'g stemming from the CT is given by 



w 



1,2,3 _ 
4,5,6 — 



E 



6 



rpp+q- 



rpj+m-r, 



rpi+j-p,p rpq+r—m,m 



i,j,m=l p,g,r=4 



?+m— r.r 

■J 

3,™ 



q-\-r—m,m 
q.r 



(27) 



where we have use the notation of w 3 ' 4 = lj% + ui^ — ^3 — <^>4- Note that analogous to the 
symmetries of the four- wave interaction coefficient Tg 1 '^ , we must similarly impose the 

7 1,2,3 
4,5,6 



following symmetry conditions on W\^r to ensure the Hamiltonian is real: 



w. 



3,1,2 
4,5,6 



4,5,6 
2,3 



(28) 



Hamiltonian (26 1 represents the the original Hamiltonian system (16), but now in the 



new canonical variable Ck- The interaction Hamiltonian has now been transformed from 
having a leading non-resonant fourth-order interaction term into one with a leading res- 
onant six-wave interaction. From the formula of the new six-wave interaction coefficient 



(27), the six- wave interaction stems from the coupling of two non-resonant four- wave 
interactions connected by a virtual wave (an illustration is presented in figure [5]). 

By substituting Hamiltonian (26) into equation (20), we get an evolution equation 
for the wave action variable Ck, 



ic k = w k c k 



1 

12 



2,3,4,5,6 



k,2,3 



fk,2,3 



4,5,6 C 2 C 3 C 4 C 5C6 04,5,6' 



(29) 



where we denote the "time" derivative of Ck as Ck = dck/dz. This equation is the starting 
point for |WT| theory. This is the six- wave analogue of the Zakharov equation describing 
four- wave interactions of water surface waves [111] . 
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5.3. The Weak Nonlinearity Expansion 

The |WT strategy exploits the separation of the linear and nonlinear timescales. 
Namely, when considering a weakly nonlinear regime, (when wave amplitudes are small), 
the linear evolution time, 

T L = — , (30) 

is much smaller than the nonlinear evolution time Tnl - the characteristic time for the 
nonlinear transfer of energy between waves. This separation of timescales allows for an 
expansion in terms of a small parameter and then an average over the fast linear dynamics 
to be taken. The result is a description of the nonlinear evolution of the wave system. 
First of all we note that the leading order nonlinear effects arise from the diagonal terms 
in the sum corresponding to the nonlinear interaction. These terms correspond to the 
trivial pairings of wave numbers in equation ( 29 ) , which do not contribute to the nonlinear 



exchange of energy between wave modes, but do provide an additional contribution to 
the frequency. Consequently, before we can proceed, we must isolate this contribution 
to ensure that the weak nonlinearity expansion will be defined in a self-consistent way 
1 104] . By separating the diagonal terms from the main summation, equation ( 29 ) can be 
expressed as 



(uj k + Q k ) c k + 



1 

12 



E 

K 



W 



k,2,3 * * 



2,3 
4,5,6 



4.5.6 



c 2 c 3 c 4 c 5 c 6 <V 5 ' 6 , 



(31) 



where 



^=oE Wk.7.8 M 2 M 2 . (32) 



2^ 

7,8 



The summation in (|31j) is now taken over the set K./ 5 6 = {k 2 , k 3 , k,j, k 5 , k 6 : k, k 2 , k 3 ^ 
k4, ks, kg}. We see that the nonlinear frequency flk acts as a correction to the linear wave 
frequency, so that the system evolves with an effective wave frequency of ujk — ujk + ^k- 
Let us define a scaled interaction representation variable 6k, so that it evolves on the 
slow nonlinear scale. This is achieved by compensating for the fast oscillating factor by 



Ok = — e 
e 



ioj k z+i J* Vt k dz' 



(33) 



Note that we have incorporated the nonlinear frequency correction f2fc into the oscillating 
factor, so that the variable 6k varies only at the timescale of the energy transfer between 
the modes. We have also divided the right-hand side by a formal small parameter e <C 1 
for easier power counting of the nonlinearity orders (thus formally defining 6k to be 
0(1)). 



Substitution of formula (33 1 into equation (31) gives 



k 12 



E 



W4 k 5 2 6 6^646566 



rk,2.3 

d 4,5,6 e 



fe,2.3 1 - 4 rz 



(34) 



■Oo + Os-Oa-Qk — fifi. As 



where w 4 'g '| = cj^ + cj 2 + w 3 — w 4 — w 5 — a; 6 and flj'gg = fife -1- 1 1.2 f i 23 — s «4 — "5 — 1 4 6- 
the nonlinear frequency correction is time-dependent, extra care should be taken when 



-u>6 and "4 5 6 



considering its contribution to (34) 
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To separate the slow and the fast timescales, let us chose an auxiliary intermediate 
time T, 

T L <T<r NL . (35) 

For small wave amplitudes, |ck| ~ e « 1, the nonlinear timescale for a six-wave process 
is of the order Tnl ~ 27r/e 8 u;fc, which will be verified a posteriori. Thus, to be specific, 
let us take T = 2w/c 4 wi ( . 

Let us seek a solution for 6 k , at the intermediate time T in the form of an e-expansion 



&k(T) 



-1.(0) 



(36) 



and then solve each e-order of bk(T) by an iterative method using evolution equation 
(34). We substitute (36) into (34), and equate e-orders. The leading O(e ) contribution 
gives 

^V)=fo k 0) (0) = 6 k (0). (37) 



Result (37), implies that at leading order, the wave amplitude bk(T) is time-independent. 
In the following order, O (e 4 ) , we get 



k,2,3 ,(0)*,(0)*,(0),(0),(0) <-k,2,3 



12 



7 1 

4,5,6 °2 



b 4 b 5 b 6 d 4,5,6 A T 



(^l + ^^ta) ■ (38) 



Here we have used that fr k (0) = and have defined 



z dz 



(39) 



and ^nlfl = n k } + ^ + - 



n 



Co) 



fif where 



yE<> 7 ; 8 8 i4 0) i 2 i4 0) i 2 - 

7,8 



(40) 



Note that, even though ftf = O (e 4 ), the factor e^'^.e J i s because T ~ e 

Therefore, the e-orders of the of the left-hand side and the right-hand side are the same. 
In other words, it would be wrong to expand the exponential function in small e here. 
The next order correction to Qk is indeed small enough (~ e 8 ) so that the exponential 
could be expanded, but its contribution would be nullified by the phase averaging (see 
below) . 

I L \ 1 (2) 

Finally, we use formulae (37) and (38) to determine & k (T): 



(0) o fc,2,3 r 



1.(2) 



1 

12 



£ 

jl-—2 , 3 if T , S 
^4,5,6' ^9, 10, 11 



2W; 



c,2,3 ( 
,5,6 ^ 



_ rk,2,3 <-2 
X d 4,5,6 °9 



2,7,8 
^9,10,11 



2,7,8 



. (0). , (0)* , (0)* , (0)* , (0) , (0) , (0) , (0) , (0) 



J 9 



(°>,7: 



fc,2,3 
•'4,5,6 



°10 °11 °4 °5 °6 °7 °8 



oyyk,2,3 W 4,7,8 



(41) 



,(0)*,(0)*,(0)*,(0)*,(0),(0),(0),(0),(0) A k,2,3 r4,7,8 p 

xb 2 & 3 b 7 b s b 5 b 6 b g b w b 11 d i56 9A0A1 ti T 



(0),7, 



(o),r, 



fc,2.3 



w 9,10,ll' w 4,5,6 



20 



where we have defined Et(x, y) as 



E T (x,y)= / A z ( x yKy dz 'dz, 



(42) 



j (o)~k.2,3 ~fc,2.3 

and ^ ; w 4 ; 5 , 6 = Wi Vs 



4,5,6 



(0) o fe,2,3 
"4.5.6- 



The symmetry of At(x) implies that Ex{x*,y) 



Et(—x, y)- With equations (37), (38) and (41 ), we have calculated the e-expansion of 6k 



up to 0(e ). This order is sufficient for the analysis contained within this review. The 
next objective is to substitute the above expansion into the definition of |GF[ thereby 
finding the solution for |GF| at the intermediate time T. 



5.4- Equation for the Generating Functional 



Recall that the one- mode amplitude |GF| (|2T|) is defined as 

2 k (A k ) 



= / P (^)^|bkl 2 



(43) 



The aim is to calculate the e-expansion of Z^ at the intermediate time, T, and sub- 
sequently, derive an evolution equation, using the results of the previous Subsection. 
Assume that we can represent the|GF|in powers of e as 



(0) 



e^> 



8^(2) 



e s Z, 



(44) 



Then to calculate the GF (43), we must average over random phases, (-)^, and random 
amplitudes, (•) j, in an initial RPA field. Once we have performed these, we can derive the 



evolution equation for Z^, and thus any of the one-mode statistical objects we require. 
The technical details of the phase and amplitude averaging are located in Appendix | B| 
After such an averaging, and after taking the large-box limit L —> oo followed by the limit 
T — > oo (corresponding to small nonlinearity e — > 0) we arrive at an evolution equation 
for the EH 



Zy = A k 77 k Z k + (A k 77 k - A k 7k) 



dZ v 
5A k 



where 



% = 



IW, 



and 



7k 



4,5,6 1 2 ^4,5,6 n 2 n 3 n 4 n 5 n 6 dk 2 dk 3 dk 4 dk 5 d'k 6 , 



i 1A ,k,2,3|2 rk,2,3 r/ k,2.3\ w i \ 

(; ; l w 4,5,6l 04,5,6 5 W,5,6) L( n 2 + "3) n 4 n 5 n 6 



(45) 



(46a) 



—77,2773 (714715 + 77477-6 + 77577.6)] (ik2dk3(ik4(ik50?k6. (46b) 
Here, we have defined the wave action spectrum 77k = {L/2tt) (Jk) with Jk = |6k| 2 - 
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5. 5. Moments of the Wave Intensity Jk 

Moments of J k , denned as M^ p) = (j£), can be easily found by differentiating 

the|GF|p times with respect to Ak and evaluating at Ak = 0: 



2 k 



(47) 



A k =0 



(p) 

Therefore, the evolution equation for can be found directly from (45 1 using formula 
p7|, giving 

(48) 



= -P7kMk P) 



2 nAP- 1 ) 



In particular, for p — 1 we have the KE for the wave action spectrum. We will consider 
the IKEI in more detail later. 



(49) 



For Gaussian wave fields we find that the p-th order moment is given by 

M L p) = pK- 



However, these solutions are invalid at small Ak and at high p, as large amplitudes, 
for which the nonlinearity is no longer weak, have a strong contribution to their value. 
Therefore, it is more suitable to consider quantities that are local in Jk, such as the 



amplitude PDF 



5.6. The One-Mode Amplitude Probability Density Function 

Recall, that the one-mode amplitude |PDF| 'Pk, can be found by applying the inverse 
Laplace transform, (22), to the GF Z^. By applying the inverse Laplace transform to 



the evolution equation for the GF (45 1, we arrive at the evolution equation for the PDF 



Vv = - 



<9s k ' 



where we have introduced a probability space flux J-k defined as 

VP* 
8sl 



Jk = - SkTkT'k + s k ?7k 



(50a) 



(50b) 



The simplest steady state solution corresponds to the zero flux scenario, J-"k = 0. This 
solution is given by 

1 _£k 

■Phom = — e n k, (51) 



where corresponds to any stationary state of thc |KE| Solution "Phomi is the Rayleigh 
distribution. Subscript horn refers to the fact that this is the solution to the homogeneous 
part of more general equation for a steady state with a constant non-zero flux, J-k — J- ^= 
0. The general solution in this case is |105j 



Pk = V ho 



part ? 



where V par t is the particular solution to equation (501, given by 



Ei 
22 



(52a) 



(52b) 



where the function Ei(x) is the exponential integral function 



Ei(x) 



-dz. 



(53) 



The particular solution is a correction due to the presence of a non-zero flux. In the 
tail, where Sk 3> ^ki we can expand "P par t m powers of n^/s^: 



region of the 



PDF 



'part — 



T 



Sk7k (s k 7k) z 



(54) 



Thus, at leading order, the PDF tail has algebraic decay ~ 0(1/ 8k). This behavior 



corresponds to the presence of strong intermittency of WT [105] . However, it was also 
noted in |105j . that if the weak nonlinearity assumption was valid, uniformly on the 
whole of Sk-space, up to Sk = oo, then the probability flux, J-^, is forced to equal zero, 
to ensure positivity of Vk and convergence of its normalization: 



Vk ds k = 1. 



(55) 



Hence, for a finite flux solution to exist, there must be a cut-off in Sk-space, for which 
the [PDF] vanishes. It arises at high wave amplitudes Sk for which the weak |WT] descrip- 
tion fails. This can be viewed as a wave breaking process, which does not allow wave 
amplitudes to exceed a critical value |105j . The wave breaking process is essential for 
the presence of WT intermittency. The sign of the flux J 7 , determines whether the tail 



of the amplitude PDF is greater (enhanced) or less (depleted) than the exponential tail 
of the Rayleigh distribution (51). From equation (54), we observe that a negative F 
implies an enhanced probability of high intensity events. Subsequently, a positive flux, 
would imply that there less probability in observing high intensity structures than 
what is expected by a Gaussian wave field. In wave turbulent systems, we expect to 
observe both kinds of behavior each realized in its own part of the k-space forming a 
|WTLC[ which will be discussed later. 

5. 7. The Kinetic Equation 

The KE is one of the key objects in WT] theory, it describes the evolution of the wave 
action k-space density, n^. The evolution of rik is governed by the equation for the first 



moment which, written in its full form using the definitions for rjk, (46a), and 7k 



( 46b ) , is expressed as 



«k 



|1«,k,2.3|2 rk.2,3 r/ fc.2,3\ 

l w 4,5,el d 4,5,6 <H W 4,5,6) n k n 2 n 3 n4n 5 n 6 



fc.2,3\ 



_1 +n, 1 -n, 1 



1 — rig 1 ) <ik2<ik3(fk4<ik5<ike 



(56) 



The integral on the right hand side of the KE (56), is known as the collision integral. 



Stationary solutions of the [KE] are solutions that make the collision integral zero. There 
exist two types of stationary solutions to the |KE| The first type are referred to as the 
thermodynamic equilibrium solutions. The thermodynamic solutions correspond to an 
equilibrated system and thus refer to an absence of a k-space flux for the conserved 
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quantities, (in our case, linear energy, 7£ 2 , and total wave action, AT). The second type 
are known as the |KZ| solutions. They correspond to non-equilibrium stationary states 
determined by the transfer of a constant non-zero k-space flux. They arise when the 
system is in the presence of forcing (source) and dissipation (sink), where there exists 
some intermediate range of scales, known as the inertial range, where neither forcing 
of dissipation influences the transfer of the cascading invariant. The discovery of the 
KZ| solutions for the |KE| has been one of the major achievements of |WT| theory, and as 
such these solutions have received a large amount of attention within the community. In 



systems that possess more than one invariant, the KZ solutions describe the transfer of 
invariants to distinct regions of k-space [3]. For many systems, these regions are usually 
the low and high wave number areas of k-space, however this is not necessarily the case 
for anisotropic wave systems |112j . The directions in which the invariants cascade can 
be discovered by considering a Fj0rtoft argument. 

5.8. The Fj0rtoft Argument 

As |fD|OWT| has two invariants, there are two [KZ] solutions of the |KE| each defined by 



a constant flux transfer of either invariant. This is analogous to 2D turbulence, where the 
enstrophy, (the integrated squared vorticity), cascades towards small scales and energy 
towards large scales (BTJ [62] . When a non-equilibrium statistical steady state is achieved 
in a weak nonlinear regime, the total energy (which is conserved) is dominated by the 
linear energy (Ji ~ ^2)- Hence, we can make the assumption that the linear energy 
is almost conserved. This is important as the linear energy is a quadratic quantity in 
ip(x, z) and allows for the application of the Fj0rtoft argument |113j . This argument was 
originally derived in the context of |2D| turbulence, and does not require any assumptions 
on the locality of wave interactions. To begin, let us define the (linear) energy flux 
P(k, t) = Pk and wave action flux Q(k,t) — Q^: 

< 57b > 

where the energy density in Fourier space is defined as £k = ^fc^k, such that H2 = 
J eje dk. Below, we will outline the Fj0rtoft argument in the context of the six-wave 
|OWT| system. 

We should assume that the system has reached a non-equilibrium statistical steady 
state, therefore the total amount of energy flux, Pk, and wave action flux, Qk, contained 
within the system must be zero, i.e. J Pk dk = and J Qk dk = respectively - this 
corresponds to the flux input equaling the flux output. Then, let the system be forced 
at a specific intermediate scale, say ky, with both energy and wave action fluxes being 
generated into the system at rates P/ and Qf. Moreover, let there exist two sinks, one 
towards small scales, say at k + ^> kf, with energy and wave action being dissipated at 
rates P+ and Q+, and one at the large scales, say at fc_ <C kf, dissipated at rates P_ and 
Q-. Further assume that in between the forcing and dissipation, there exist two distinct 
inertial ranges where neither forcing or dissipation takes effect. In the weakly nonlinear 
regime, the energy flux is related to the wave action flux by Pk ~ UkQk — k 2 Q^. In a 
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non-equilibrium statistical steady state, the energy and wave action balance implies that 



Pf = P- 

Qf = Q- 

and therefore, we approximately have 



P 



f 



P_ 
P+ 



-P+, 

-Q+, 

k}Q f , 

k 2 _Q-, 
k 2 + Q+. 



Subsequently, the balance equations ( 58 ) imply 

',.2 
V 



k 2 Q f 



Q f 



k 2 _Q_+k 2 + Q+, 
Q- + Q+. 



(58a) 
(58b) 



(59a) 
(59b) 
(59c) 

(60a) 
(60b) 



Re-arranging equations ( 60 ) enables us to predict at which rates the energy and wave 



action fluxes are dissipated at the two sinks. From equations (60 1 we obtain 

k 2 



Q- 



k)_ 

k 2 - k 2 Qh 



k 2 



ft 
k 



2 Qf- 



(61a) 
(61b) 



If we consider the region around large scales, fc_ <C kt < fc+, then equation (61a) implies 

k 2 Q f « k 2 + Q+, (62) 

Furthermore, considering 



i.e. that energy is mostly absorbed at the region around k 
the region around small scales, k_ < kj <C k + , (61b) implies that 



(63) 



i.e. that wave action is mostly absorbed at regions around fc_. Ultimately, if we force the 
system at an intermediate scale, where there exists two inertial ranges either side of kf, 
we should have the majority of the energy flowing towards small scales and the majority 
of the wave action flowing towards large scales. This determines the dual cascade picture 
of the six- wave system illustrated in figure [6| 

5.9. The Zakharov Transform and the Kinetic Equation Solutions 

To formally derive the thermodynamic and |KZ| solutions of the |KE| we will use the 
Zakharov transform (ZT). This requires that the interaction coefficients of the system 



are scale invariant. Scale invariance of an interaction coefficient is reflected by its self- 
similar form when the length scales are multiplied by a common factor, i.e. for any real 
number A, we say that an interaction coefficient is scale invariant with a homogeneity 
coefficient f) € R, if 



W(Aki, Ak 2 , Ak 3 , Ak 4 , Ak 5 , Ak 6 ) = A^W(ki, k 2 , k 3 , k 4 , k 5 ,k 6 ). 
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(64) 




Dissipation 



Dissipation 



k_ k f k + 

Figure 6: A Graphic to show the development of the the dual cascade regime for |lD|OWT| 



Moreover, the frequency ujk must also possess the scale invariant property, i.e. 



w(Ak) = A a u;(k), 



(65) 



with some a 6 
we see that a = 



For 



OWT this is indeed the case. Since for OWT we have w(k) cx k 



2. Let us now seek solutions of the |KE| with a power-law form 
Uk = CkT x \ 



(66) 



where C is the constant prefactor of the spectrum, whose physical dimension is deter- 
mined by the dimensional quantities within the system, and x is the exponent of the 
spectrum. 

An informal way in determining the exponent x of the |KZ| and thermodynamic solu- 
tions is to apply a dimensional analysis argument. For the thermodynamic equilibrium 
solutions, we assume a zero flux, i.e. that both and are scale independent. Con- 



versely, for the derivation of the KZ solutions, we want to consider a wave action density 
scaling that implies a constant flux of the cascading invariant. This is achieved by con- 
sidering Pj., Qk cx k in equations 



Pk 

<2k 



de k , 

dz 



oz 



(67a) 
(67b) 



derived from equations ( 57 1 , using the power-law ansatz ( 66 1 and the KE ( 56 ) . However 



this method does not allow for the evaluation of the prefactor of the spectrum in ( 66 1 



Therefore, we will now describe the formal way of calculating the |KE| solutions by using 
the |ZT| The |ZT| expresses the |KE| in such a way that it overlaps sub-regions of the |KE| s 
domain of integration, thus at each solution, the integrand of the collision integral is set 
exactly to zero over the whole domain [3], The ZT takes advantage of the symmetries 



possessed within the |KE| by a change of variables. In our case, this results in splitting 
the domain of the |KE| into six sub-regions. 

Applicability of the |ZT| requires the locality of wave interactions, namely that only 
waves with a similar magnitude of wave number will interact. The criterion of locality is 
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equivalent to the convergence of the collision integral. Locality of these solutions will be 
checked in the following Section. 

The |ZT| is a change of variables on specific sub- regions of the domain, one such sub- 
region is transformed by 



k 2 

k 2 = ~~, k 3 
k 2 



kkj 

k 2 : 



kk 4 



k 5 = -= — and k 6 
k 2 



kkg 

k 2 : 



(68) 



with the Jacobian of the transformation J = — ^k/k 2 J . We must apply four simi- 
lar transformations, to each of the remaining sub-regions (these are presented in Ap- 
pendix C ) . 



Using the scale invariant properties of the interaction coefficients and frequency, and 
the fact that a Dirac delta function scales as 



S((Xk) a ) = \- a 5(k a ), 



(69) 



the |ZT| implies that the KE can be expressed as 

J |W 4 k ; 5 2 ; 6 3 | 2 |kk 2 k 3 k 4 k 5 k 6 p {k x + kl + k* -k%- kl - k%} 



C 5 e 8 7 



n k = 



k^ y 



k^ y 



k 4 ^ y 



k 



^4,5,6 ^(^4,5,6) rfk 2 dk 3 dk 4 dk 5 dk 6 , 



(70) 



where we have omitted the tildes and y — 5x — 3 — 2/3. We see that if x = or x = 2, 
the integrand will vanish as we have zero by cancellation in either the curly brackets or 
in the square ones (taking into account the Dirac delta function involving frequencies). 
In particular, when x = or x = 2, the solutions will correspond to thermodynamic 
equilibria of the equipartition of the energy and the wave action respectively: 



Tij^ — Cu k 



l k — Cjfk , 



(71a) 
(71b) 



Solutions ( 71 1 correspond to zero flux states - in fact both energy and wave action 
fluxes, Pk and Q k are identically equal to zero on both equilibrium solutions. Spectra 



( 71 ) are two limiting cases, in low and high wave number regions, of the more general 



thermodynamic equilibrium Rayleigh- Jeans solution n^ J , where 



,ru 



LU k + fj, 



(72) 



Here T c is the characteristic temperature of the system, and n is a chemical potential. 
In addition to these thermodynamic solutions, our system possesses two non-equilibrium 



KZ solutions. The KZ solutions are obtained from equation (70) when either y = or 
y = 2. When cither condition is met, the integrand in equation (70) vanishes due to can- 



cellation in the square bracket. The corresponding solution for the direct energy cascade 
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from low to high wave numbers is obtained when y 
action spectrum scaling: 



= C n k 



"H-2;-i 



2, this gives the following wave 

(73a) 



The wave action spectrum ( 73a ) implies that the energy flux Pk is k- independent and 



non-zero. The solution for the inverse wave action cascade from high to low wave numbers 
is obtained when y = 0, and is of the form 



(73b) 



On each KZ solution, the respective flux is a non-zero constant - reflecting the Kol- 
mogorov scenario, whilst the flux of the other invariant is absent. However, we emphasize 
that the |KZ| solutions are only valid if they correspond to local wave interactions - an 
assumption of the |ZT| 

The main contribution to the six- wave dynamics in the |LWE[ after application of the 



CT and expansion in small kl^ (see Appendix D), is k-independent: 



1,2,3 
4,5,6 



(74) 



This implies that the homogeneity coefficient for the LWE is (3 = 



Therefore, equations ( 73 1 imply that the KZ solutions for the LWE are given by 

4/5 



= a 



n 



K' 5 k- 



(75a) 



n k 



C L 



M 



4/5 



Qi /5 fc- 3/5 , 



(75b) 



where L n^, is the KZ spectrum for the direct energy cascade, and L n^ is the inverse 



wave action spectrum. Here, and Ch are dimensionless constant pre-factors of the 
spectra. 



Conversely, for the 



SWE 



we see that the four-wave interaction coefficient, S T^'^, is 



scale invaria nt and scales as oc k 2 . Formula (27), implies that the homogeneity coeffi- 
cient for the SWE is (3 s — —6. We calculated the explicit expression using Mathematica, 

4,5,6 



and confirmed that /3 



s _ 



-6. However we must omit the explicit expression for W, 



in this review as it is extremely long. 

Therefore, the SWE] yields the following KZ solutions: 



S 7i 



III, 




4/5 



P k 1/5 fc 7 /5, 



4/5 



qT fc 9/5 



(76a) 



(76b) 



where, and Cfj- are dimensionless constants. 



Spectra ( 75 ) and ( 76 1 are valid solutions only if they correspond to local wave inter- 



actions. Therefore we must check that the collision integral converges on these spectra. 
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5. 10. Locality of the Kolmogorov-Zakharov Solutions 

The |KZ| solutions can only be realized if they correspond to local, in k-space, wave 
interactions. This entails checking that the collision integral converges, when the wave 



action density is of KZ type, (73). Our strategy for determining the locality of the KZ 



solutions is to check the convergence of the collision integral when one wave number 
vanishes, or when one wave number diverges to ±oo. These limits correspond to the 



convergence of the collision integral in the infrared ( IR I and ultraviolet ( UV ) regions of 
k-space. Although the collision integral is five-dimensional (5D), the two Dirac delta 



functions, relating to the six-wave resonance condition, imply that integration is over 
a 3D surface within the 5D domain. For linear frequencies of the form, ujf. oc k 2 , we 



can parametrize the six- wave resonance condition (A.7|, which subsequently allows us to 



neglect two of the integrations in ( 56 1 



The details of the locality analysis is situated in Appendix | E| We performed the 
analysis using the Mathematica package, which allowed us to handle and simplify the vast 
number terms resulting form the |CT| Assuming that the six- wave interaction coefficient 
has the following scaling as — > 0, 



lim W 



k,l,2 
4,5,6 



k 6 
K 6l 



(77) 



where ^eR. Then we find that we have |IR| convergence of the collision integral if 

x < 1 + 2£, (78) 



is satisfied. Similarly, by assuming that the six-wave interaction coefficient scales with 

lim w 4 k ; 5 2 ; 6 3 cx kl (79) 



/ug as kt 



where r/ G R, we find that the criterion for |UV| convergence of the collision integral is 

7] < x. (80) 

We now ch eck fo r convergence in the |OWT| models. We begin by investigating the 
The six-wave interaction coefficient L yv]'5 g was shown (at leading 



LWE 



locality on the 

ord er) to be constant, with the constant given in equation (74). This im plie s that in the 



IR 



limit L Wlfe remains constant, i.e. £ = 0. Hence, the condition for 



Becomes 



IR 



x < 1. 



convergence 



(81) 



Due to relation (74| being constant, we also have that r\ = 0, and therefore, the conver- 
gence condition for 



UV| region is given by 



< x. 



(82) 



Therefore, the |LWE| has local |KZ| spectra if their exponents are within the following 
region: 

< x < 1. (83) 



For the LWE s KZ solutions ( 75 ), we have that the direct cascade of energy has the expo- 
nent x = 1, which implies divergence. However, because this exponent is at the boundary 
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of the convergence region, we have a slow logarithmic divergence of the collision integral. 



This implies that by making a logarithmic correction to the wave action spectrum (75a), 
we can prevent divergence of the collision integral - we will consider this in the next 
Subsection. The inverse cascade of wave action, with x = 3/5, implies convergence of 
the collision integral. 

Consideration of the |SWE| in the [TO] region, by appropriate Taylor expansion around 
a vanishing kg using Mathematica, reveals that the six-wave interaction coefficient, 
s W45'g, behaves as 



15m S <; 6 3 oc fca 1 , 

ke-s-0 



(84) 



giving £ = — 1. Therefore, the condition for |IR| convergence of the [KZ| solutions becomes 



x < -1. 



(85) 



Similarly, using Mathematica and expanding the interaction coefficient >y^§§ 
limit where kg — > oo, we find that 



s Wh 2 ' R 3 in the 



lim s W^oc*§, 



(86) 



thus r\ = 0. This implies that the UV condition for convergence is the same as that for 
the LWE relation (82). Therefore, there does not exist any region of convergence for 



the collision integral of the SWE To be specific, both KZ solutions of the SWE (76), 
where x = —7/5 and x = —9/5 for the direct and inverse spectra respectively, we have IR 
convergence and |UV | di vergence. Therefore, both short-wave KZ spectra are non-local. 
Non-locality of the KZ solutions implies that the local wave interaction assumption is 
incorrect, and thus the approach taken to predict these spectra is invalid. However, the 
development of a non-local theory for the |SWE| may yield further insight. 



5.11. Logarithmic Correction to the Direct Energy Spectrum 

In the previous Subsection, the direct energy cascade in the LWE 
to be marginally divergent in the |IR| limit. This is to say, that the collision integral 
diverges at a logarithmic rate in the limit of one vanishing wave number. However, by 



(75a), was shown 



introducing a logarithmic dependence to the KZ solution, we can produce a convergent 
collision integral. Following Kraichnan's argument 6 1 62 . let us assume a correction of 
the form: 

/ \ 4 / 5 

i>k\ „i/6,__i,_-„ (w)j (87) 



C L 



hf 



where y is some constant to be found and £ is the scale of at which energy is injected. 
The exponent of the logarithmic power law, y, is calculated by assuming that the energy 
flux, Pk, remains k-indcpcndcnt. Subsequently, the energy flux can be expressed as 



Pu = 



at 



dk: 



k- L \n- 5y (k£) dk, 
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where we have taken into account that the collision integral ( 56 1 , with interaction coef- 
ficients (18 1, scales as oc k 2 n^. Therefore, remains k-independenl[^]when y = 1/5. 
This implies that the logarithmically corrected direct energy [KZ] spectrum is given by 



4/5 



hi 



-1/5 



(M). 



(89) 



Spectrum (89) now produces a convergent collision integral for the LWE and is subse 



qucntly a valid |KZ| solution. 

5.12. Linear and Nonlinear Transfer Times and The Critical Balance Regime 

In this Section, we estimate the nonlinear transfer times of the |KZ| cascades and 



discuss the critical balance ( CB I phenomenon of strong WT The separation of timescales 
is essential for the construction of the WT closure (see Subsection 5.3 1. The linear 
timescale was previously defined in equation ( 30 ) . It can be noted that from equation 
(35) the nonlinear timescale, Tnl, is 



C(l/e 8 ). From the KE we can define Tnl as 



dt 



(90) 



In general, this expression is /c-dependent. |WT| theory is applicable when there is a large 
separation between the linear and nonlinear timescales, i.e. 



Tnl 



< 1. 



(91) 



The estimation of Tnl, with respect to k, can be achieved using the KE (56), giving 

T NL oc fc 4 *- 2 ^ 2 . 



Therefore, the condition of applicability (91 ) can be written as 



T L 
Tnl 



k 



-4x+2/3 



< 1. 



(92) 



(93) 



Relation (|93|) can be violated in either of the limits k 
sign of 



or k — > oo, depending on the 



1, then the KE approach breaks down and we are in 



-Ax + 2,3. When T L /T NL 
a strong |WT] regime. 

To describe strong [WT| in some physical systems the concept of |CB| was suggested, 
first in magneto-hydrodynamics (MHD) turbulence |32j and then later in rotating and 
stratified turbulence |114) and BEC 47J. CB is defined as a turbulent state where the 
nonlinear evolution time Tnl is of the same order as the linear wave period Tl over a 
large range of scales . A |CB| scaling for the wave action density, n^, can be made by 
equating Tstl with Tl, on a scale by scale basis, i.e. 



^ oc k- 4x+2 ? ~ 1, 
Tstl 



(94) 



3 The energy flux is actually proportional to oc ln(ln(fc^)). 
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implying a |CB| spectrum of 



C C Bk 2 . 



(95) 



where Ccb is a constant prefactor of the spectrum. 

For |1D||0WT| we est ima te the ratio b etween the linear and nonlinear timescales 
using the 



KE 



of the 



LWE 



and 



SWE 



0, 



For thelLWEl with homogeneity coefficient (3 L 
implies that for both the direct (x = 1) and inverse (x — 3/5) KZ spectra, the WT 
criterion (91) will gradually be violated as k — > 0. For the direct cascade, the largest 



scale is at the forcing. Thus, if |WT| is weak at the forcing scale, it will remain weak 
everywhere along the cascade. In the inverse cascade, even for weak forcing, the strength 
of |WT| is increasing along the cascade toward the lowest wave numbers, resulting in the 
break down of WT| theory and the formation of nonline ar co herent structures. For the 



SWE with homogeneity coefficient f3 



s _ 



6, and with KZ exponents x = —7/5 and 



9/5 for the direct and inverse KZ spectra respectively, we have a similar situation 



where condition (91 ) gets violated as the inverse cascade progress to small wave numbers. 



5.13. The Differential Approximation for the Kinetic Equation 



The differential approximation model ( DAM I is an approximation of the KE by as- 
suming strongly local wave interactions. This enables the construction of a partial differ- 
ential equation for the evolution of the wave action density, n^. We stress that the DAM 



is an approximation and thus can only be justified when the |KZ| solutions are proved local, 
i.e. the |DAM| is only applicable to the|KE]of thc |LWE| Therefore, we find that the |DAM| 
is only applicable to the KE| of the |LWE The |DAM contains both the thermodynamic 
and non-equilibrium |KZ solutions of the |KE| and can be further simplified to reduced 
DAM j, which only consider a subset of these solutions. The usefulness of the DAM can be 



shown with derivation of exact analytical solutions, expressions for the fluxes and finally 
by its computational simplicity. The |DAM| can also be adapted to classical turbulence 
theory, where it is known as the Leith model [1 151 11161 1117] . This has led to the DAM be- 
ing used extensively in |WT|a nd classical turbulence [XX51 IXX91 11201 1X211 1X221 HX1 1X251 IXX0] . 
The derivation of the |DAM| for the |LWE| can be found in Appendix [T] The |DAM| for 
l~D|OWT|is given by 



duj 2 



(96) 



This equation contains the same solution (thermodynamic and non-equilibrium) as the 



KE The thermodynamic solution, n u — T c /(uj + /i) is found when 



d 2 



dui 2 V n, 



1 



(97) 



is identically equal to zero. The energy and wave action fluxes, and Q U1 can be 
derived from R M by the following formulae: 



P* 



R^ -uj 
dRy 
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du ' 



(98a) 
(98b) 



Consequently, both P u and Q M vanis h up on the thermodynamic solution. By assuming 



the same power-law scaling as for the |KE| (|66| , n w = Clo x / 2 , we can calculate formulae 
for the behavior of the |DAM] with respect to to and the exponent x: 



R{uj,x,y) 



= S a C 5 
= S C 5 



(f-l) »(!/-!)"-»-* 
Q(w,x,y) = S C 5 |(|-l)(y-l)^, 



P(w,x,y) = S Q C ; 



X /X 

2 

a; f .< 

2 
.r 



-v+i 



-y+i 



where y — - 
x = 0, x = 2, 



I. Stationary solutions for the DAM are seen when n u 
I) 



(99a) 
(99b) 
(99c) 
(99d) 
0, i.e. when 



or y = 1, yielding the thermodynamic solutions (71) and the 



KZ solutions of the LWE ( 75 1 . Relations ( 99c ) and ( 99d I enable for the calculation of 



the sign of the KZ fluxes. First, we observe that both fluxes vanish upon reaching the 
thermodynamic solutions, when x = and x — 2, and second, that on the KZ solutions, 
the flux of the non-cascading invariant is zero. 



The Fj0rtoft argument of Subsection |5.8| showed the directions of k-space in which 
energy and wave action are permitted to flow. However, the direction of the invariant 



cascade in the KZ solution may sometimes contradict Fj0rtoft's argument in which case 
they cannot be matched to any physical forcing and dissipation, and hence, they are not 
realizable. Therefore, it is essential that the direction of the flux agrees with the analysis 
of Fj0rtoft's argument. 

We use the formulation of the |D AM] in Subsection [5T3 to determine the directions of 



the energy and wave action fluxes Pk and Qk within the LWE By plotting formulae ( 99c 



and (99d) for Pk and Qk against the exponent of the KZ solution, x, we can determine 



the sign of the fluxes. When the exponent is that of both thermodynamic solutions (71 1 
then both fluxes should vanish. 



In figure [7J we plot equations (99), for P k and Qk- we observe that at x — 0, 2 
corresponding to the two thermodynamic equilibrium solutions, we have both the energy 
and wave action fluxes identically zero. At x = 3/5, corresponding to the wave action 
cascade, we have that the energy flux Pk is zero, whilst the wave action flux Qk is 
positive. When x = 1, the exponent for the energy cascade, we have the wave action 
flux Qk is zero, and the energy flux Pk negative. However, to agree with Fj0rtoft's 
argument, when x = 1, the energy flux must be positive and the wave action flux zero. 



Similarly, when the exponent agrees with the KZ exponent of the inverse cascade, (75b) 



Fj0rtoft's argument implies that the wave action flux must be negative. Thus, we should 
not expect realizability of the pure KZ solutions. The authors of 48 suggest that in 



such situations, we might expect to observe a finite temperature cascade solution. This 



solution is predominantly a thermodynamic solution similar to ( 72 ) , but with a finite 



flux contribution P(wfe). This results in a mixed solution of thermodynamic equilibrium 
and non-equilibrium parts of the form 



T, 



n k = 



0J k + U + P(Wfc) ' 



(100) 



with F(uj k ) <C oj k , (i. 
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5.14- Modulational Instability and Solitons in the Long-Wave Equation 

The similarity of the LWE| Q to the integrable ID focusing NLSE| ^ means that we 
should expect not only random waves but also soliton-like coherent structures. During the 
inverse cascade, solitons appear naturally as wave action reaches the low wave numbers. 
Indeed, the KE description (56) becomes invalid, and the MI of the wave field occurs 



resulting in the filamentation of light and its condensation into solitons. 

Solitons are permanent localized (in physical space) structures resulting from the 
balance between the dispersion of the wave and the nonlinear wave breaking effect of 
the focusing medium. Solitons appear in several nonlinear wave equations, with the 
ID NLSE| one example. The ID NLSE can be solved exactly by the inverse scattering 
transform. Analytical soliton solutions in this case have the form 



ipsoi (x, z) = 2 A ip sech ( Ax 



ABz 



C\ exp 



(A 2 - B 2 ) z 

Bx + ± >— + C 2 

2q 



(101) 



where A, B, C\ and C 2 are constants. 

Unlike equation the |LWE[ 0, is non-integrable, and thus, will not possess exact 
analytical solutions of form ( jlOl] ). However, the LWE s deviation from NLSE is small, 
so we would expect to observe coherent structures of the LWE to be close in shape to 
(101). Moreover, solitons of integrable systems possess the property of passing through 



one another without changing shape. However, when integrability is slightly violated, we 
expect to observe weak nonlinear interactions between solitons. This may be observed as 
possible merging events, oscillations in profile or interactions with weak random waves. 
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Formally, [lD|NLSE| solitons are denned as solutions to ^ of the form 



Vz 



(102) 



i.e. that the structure of the soliton is preserved as it propagates along the z-axis by some 
velocity V . The Fourier transform of <fi(x, z) satisfies the following integral equation: 



1 



2(fl-kV + u} k ) 



T, 



k,2 /* 

3,4 92<P3<P4a 3 4 



(103) 



with T 3 4 , in the case of the 



ID 



given by the first term of equation (18 a), i.e 



NLSE 

constant. If the denominator ot 0k m equation ( 103 ) has no zeros for k € 
may be present in the system. 



, then solitons 



If in equation ( 103 1 the numerator and denominator 



contain only one zero each at k = then solitons may also exist. For |OWT| when 
cjfc = fc 2 , there exist soliton parameters fl and V for which the denominator is not zero. 
This corresponds to the cases when on the (k, cj)-plane, the line uj = kV — lies below 
the parabola u) = uik = k 2 without intersecting it. Such parameters correspond to non- 
radiating solitons, whereas for lines intersecting the parabola we expect structures whose 
shapes evolve due to interactions with weak waves. In one scenario, such structures 
would decay via radiating weak waves (cf. Cherenkov radiation). In fact, we will see 
later that the dominant interactions are in the opposite direction: the waves "condense" 
into solitons with parameters fl and V such that the line u> — hV — is tangentially 
touching the parabola uj — ujk- Subsequent evolution brings these lines further down 
thus decoupling the solitons from the weak waves. 

In figure[8j we plot the profile of the NLSE soliton ( 101 ), with the corresponding wave 
action spectrum given in figure [9] We observe a scaling of oc k° towards low wave 
numbers, with a decline of the spectrum at large k. However, a nk cx k° scaling is also 



observed from the equipartition of wave action ( 71 1 . Therefore, by mere observation of 



the wave action spectrum, n^, it will be difficult to determine if the wave field is comprised 
of equilibrated random waves or if is in the presence of solitons. To distinguish between 
these two states, we can numerically produce a (k, w)-plot, that involves an additional 
Fourier transform over a time window. This method separates the random waves from 
the coherent component by observation of the dispersion relation, ujk B51 Wf\ . 

The inverse cascade of photons is an essential step in the creation of solitons from 
a weakly nonlinear random wave field. The cascade provides the means, via nonlinear 
wave interaction, for wave action to reach large scales. As the inverse cascade devel- 
ops, the nonlincarity of the system increases and a condensate forms, resulting in the 
dynamics of random waves deviating from the linear dispersion relation to one that is 



Bogoliubov modified. In thc |LWE] this condensate is unstable, and undergoes [MI] result- 
ing in the formation of solitons. To derive an expression for the Bogoliubov dispersion 
relation, we must first expand the wave function ip(x, z) around a homogeneous conden- 
sate solution. The condensate solution describes the behavior of the zeroth mode. The 
Bogoliubov dispersion relation then comprises of the wave frequency for the disturbances 
upon the condensate. The derivation of the Bogoliubov dispersion relation can be found 
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Figure 8: Intensity profile of the |lD|NLSE| soliton profile given in equation 
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Figure 9: The corresponding wave action spectrum for the soliton profile in figure [8] 



in Appendix | G[ The dispersion relation is given by 



Io_ 

2^1} 



4> 2 



^2,2 



(104) 
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where io = IV'ol 2 is the intensity of the condensate. We observe that the Bogoliubov 
dispersion relation (104) is the linear wave frequency ujf. = k 2 (when Iq — 0) with a 



correction due to the formation of the condensate. Frequency (104) becomes complex 



when the term inside the square root turns negative. This gives a condition for the MI 

io 



(ft 



(105) 



It is MI that precedes and assists the formation of solitons when WT fails to remain 
weak at large scales. Note that the large scale condensate is not necessarily uniform 
in the physical space - it may contain a range of excited large-scale modes. However, 
such a condensate will be felt as quasi-uniform by modes with k significantly exceeding 
the lowest wave numbers. Thus, it is best to generalize the definition of the condensate 
intensity. In fact, we found that the best agreement of the Bogoliubov dispersion relation 



to the numerical (k, w)-plot (figure 21 ) was achieved when we take Iq to equal the intensity 
of the soliton component, i.e. the wave action contained in the region — w max /2 < ui < 0. 
Note that another classical example of |MI| was also discovered in the context of water 
waves, where it is known as the Bcnjamin-Feir instability |124j . 

Consideration of a single 1D|NLSE soliton defined by equation ( 101 1, in an infinitely 
long domain, allows for the analytical computation of the conserved quantities of thc |lD| 
|NLSE| system. This will provide some insight into the behavior of solitons in the |LWE| 
Calculation of the two energy contributions in the |NLSE| Hamiltonian are given by 



Ho 







I —oo 


dx 2 



dx 



8A^ 2 l 2 (A 2 + 3B 2 



Aip 2 l 2 J- 



|Vw| 4 dx = - 



16A 3 ^ 2 l 2 



(106a) 
(106b) 



Although, exact |NLSE| solitons do not interact, solitons in the |LWE| will because of 
the additional non-integrable term. Therefore, we will most probably observe behavior 
such as growing, shrinking, merging and collapsing. If we consider the A-dependcncc 



of equations (106), then we observe that the nonlinear energy of the soliton increases 
at twice the rate (with respect to A) compared to the linear energy. However, for the 
system to conserve energy, the linear and nonlinear energies must remain balanced for 



all time (with their sum equal to some initial energy, E ). From (106) this initial energy 
is given by 

^ 2 l 2 A (W 2 - A 2 ) 
H sol = E = - . (107) 



For the same energy, equation (107) implies the relation 



B 



Eo 
3A 



(108) 



i.e. for the soliton to grow in amplitude A, it must also increase its speed ~ B/q, in 
order to compensate for the balance between the linear and nonlinear energies. However, 
this would increase the wave action J\f and the linear momentum M - another conserved 
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quantity of the system. Therefore, an important question to ask is, can two solitons 
merge, while wave action, energy, and linear momentum remain conserved? 

The answer to this is yes, but only if the merging of two solitons generate waves. Let 
us prove this fact. Consider two solitons with parameters {^4i,i?i} and {^42,-82}, that 
defines their amplitude and speed respectively, and then assume that the two solitons 



merge into one soliton defined by {A 3 , B 3 }, see figure 10 for an illustration, 
Before Merging ♦ After Merging 




Figure 10: Graphic to show the before and after states from a soliton merger. On the left image we 
see two solitons with parameters {Ai, B\\ and {A2, B2}, while after we see only one soliton defined by 
{A3, B3} but with the generation of waves, with the wave energy being determined by equation JTl3t . 

The wave action J\f of a soliton is defined as 



A4oi= / |Vw| 2 dx = AA^% 



and its momentum Ai is given by 



M soi = 



sol- 



sol 



dip: 



sol 



dx = SABi]} 2 l\. 



(109) 



(110) 



2 J -00 dx ^ so1 dx 

Conservation of the total waveaction and the total momentum of two solitons imply that 



^3 - 

B 3 = 



Ai + A 2 , 
A X B X + A 2 B 2 



(111a) 
(111b) 



Ai + A 2 ' 

must hold during the merging event. Now if we consider the total energy conservation, 
but include a term that describes the generation of energy associated to waves H wa , vcs (it 
is reasonable to assume that low amplitude waves do not significantly contribute to the 
total wave action and the momentum of the system), so that the conservation of energy 
implies 

'Hsoll + ^sol2 = 'HsolS waves ? 

(112) 

where H so u is the energy defined by equation (107) of soliton i. Then using (111), (112) 
and (107), we have that the energy that is generated into waves is given by 

^i> 2 l\A l A 2 \{Ai + A 2 f + (Bi - B 2 



H 



A 1 
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An 



(113) 



Equation (113) states the amount of energy that is transferred to waves from the soliton 
merging event and is solely determined by the parameters of the initial (pre-merging) 
solitons. Furthermore, equation (113) implies that there must be waves generated from 
the merging of two solitons. If there were no waves created from the interaction, then the 
only solution to (113) would imply A\ = or A2 — 0, i.e. one of the two initial solitons 
is absent. Therefore, we can expect the behavior of the |LWE| to be as follows: As wave 
action reaches the largest scale by the inverse cascade, solitons begin to emerge via |MI| 
As these solitons evolve, they grow in size by absorbing wave action from waves and by 
merging with other solitons. But by merging with other solitons, the interaction emits 
low intensity random linear waves back into the system. This process proceeds until 
there is only a single soliton remaining surrounded by a field of random linear waves. 



6. The Numerical Method 

For the numerical computation of the |LWE| and the |SWE| we implement a pseudo- 
spectral method [1251 1126] , with a resolution of N = 2 12 Fourier modes. The scheme 
utilizes the fast Fourier transform to convert physical space vectors into their Fourier 
representation, where differentiation of the wave field is transformed to algebraic multi- 
plication by wave number k. The linear terms arc solved exactly in Fourier space and then 
converted back into physical space with the aid of integrating factors, to greatly improve 
the numerical stability. Nonlinear terms arc more complicated. Multiplications involving 
tp(x, z), have to be performed in physical space. Moreover, multiplication produces alias- 
ing errors, due to the periodicity of the Fourier series expansion of the solution. These 
errors can be removed by artificially padding the Fourier mode representation with wave 
number modes of zero amplitude at either ends of k-space. The amount of de-aliasing 
is subject to the degree of nonlinearity of the equation. For us, this entails that half of 
the wave number resolution for the solution needs to be zero - a quarter at each end of 
k-space |126j . 

We evolve the simulations in time by using the fourth order Runge-Kutta method 
|126) . using a time step that is small enough so that it satisfies the Courant-Friedrichs- 



Lewy ( CFL ) condition of 

max ( ^ I At < Ax, (114) 



k V dk 

where Aa; = L/N is the spacing of our spatial grid and L = 32ir, is the length of our 
periodic box. Numerically we compute the non-dimensionalized equation given by 

i^ = -^+N(iP)+i(F k -D±), (115) 



where N(?p) represents the nonlinear part of our equations. Appendix 



H includes 



the derivation of the non-dimensional versions of both the ILWEl and ISWE1 The non 



dimensional LWE contains an adjustable parameter a, defined in Appendix | H| For the 



decaying setup, we set a — 128, so that the numerical simulation is in the same regime 
as the experiment. For the forced and dissipated cases, we set a = 1024, so that the 
long-wave limit is better realized. = F(k) is the forcing profile, where energy and 
wave action are injected into the system. We define this in Fourier space, over a specific 
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range of wave numbers. For the direct cascade simulations we use a forcing profile given 
by 

^direct = f Acxp(i6» k ) if 9 < 2& < 11 / n6 N 

k 1 otherwise, 

where A is the amplitude of forcing and #k is a random variable chosen from a uniform 
distribution on [— 7T, n) at each wave number and at each time step and Afc = 2tt/L is 
the Fourier grid spacing. For the inverse cascade simulations, we apply forcing over a 
small range of wave numbers situated in the high wave number region. However, we must 
allow for some Fourier modes with wave numbers larger than the forcing ones, because 
some direct cascade range is necessary for the development of the inverse cascade. Our 
forcing profile for the inverse cascade simulations is given by 

^inverse = f A CX P (l8 k ) if g - 10 < ^ < $ + 10 

k | otherwise. 

In all simulations, we dissipate at high wave numbers by adding a hyper-viscosity term 
oc k 4 ipk. We use the fourth power of k, so that the dissipation profile is not too steep 
for the formation of a bottleneck, or too shallow as to prevent a large enough inertial 
range developing. At low wave numbers, we use two types of dissipation profile. Firstly, 
we can use a hypo- viscosity term oc k 4 -0k- However, this type of dissipation profile in 
the |3D||NLSE| model has led to the |WT description becoming invalid and a |CB| regime 
to develop [37]. In these situations, the CB scenario is avoided by the use of low wave 
number friction. 

To summarize, our numerical dissipation profile Z?k = D (k), removes wave action 
and energy from the system at low and high wave numbers. For clarity, we split D k into 
the low and high wave number contributions and D(? . At high wave numbers our 
hyper-viscosity profile is defined as 

D" = ^hypor k 4 Vk, (H8a) 

where ^h ype r is the coefficient for the rate of dissipation. If we apply hypo-viscosity at 
low wave numbers, then is given by 



D 



L / ^hypo k 4 ?Ak if k ^ 

V-k = if k = 0, 



k 



where ^h yP o is the coefficient for the rate of dissipation at low wave numbers. However, 
in situations where we use friction, then Z?£ is defined as 

D L _ f ^friction 0k if < fc < 6 

k 1 otherwise, 

where fMction is the rate of friction dissipation. 



7. Results 

We present the experimental and numerical results for |OWT] We divide this Section 
into two main parts, firstly, we consider the system described by the |LWE[ in both ex- 
perimental and numerical aspects. We begin this Subsection by considering the decaying 
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setup with condensation at the lowest wave numbers. Then we move on to the numerical 
simulations of the |LWE| in the forced and dissipated regime. In the second Subsection, 
we present the numerical results for the |SWE| 



7.1. The Long-Wave Equation 

7.1.1. The Decaying Inverse Cascade with Condensation 

The decaying simulation leading to photon condensation was originally reported by 



us in 49 . Both experimental and numerical setups are configured for decaying OWT 



where an initial condition is set up and allowed to develop in the absence of any forcing 
or artificial dissipation. We perform the numerical simulation with the same parameters 
as the experiment, and present results in dimensional units for comparison. 

The inverse cascade spectrum is of a finite capacity type, in a sense that only a finite 
amount of the cascading invariant (wave action in this case) is needed to fill the infinite 
inertial range. Indeed, this is determined as the integral of the wave action spectrum 
converges at k = 0, i.e. 

n k dk cx / fc~ 3/5 dk < oo. (119) 



For finite capacity spectra, the turbulent systems have a long transient (on its way to 
the final thermal equilibrium state) in which the scaling is of the KZ type (provided this 
spectrum is realizable, i.e. local and agrees with Fj0rtoft). This is because the initial 
condition serves as a huge reservoir of the cascading invariant. 

Experimentally, the initial condition is setup by injecting photons at small spatial 
scales by modulating the intensity of the input beam with a patterned intensity mask. 
We randomize the phases by the use of a phase modulator. This is made by creating a 
random distribution of diffusing spots with the average size ~ 35 (xm through the |SLM| 
This is done in order to create an initial condition close to an |RPA| wave field required 
by the theory. The numerical initial condition is more idealized. We restrict the initial 
condition to a localized region at small scales. The initial profile is given by 



W j ~ 1 ^k(O) - otherwise, [LZU) 

where A = 1.573 x 10 -8 is the amplitude of forcing and 8^ is a random variable chosen 
from a uniform distribution on [— tt, it) at each wave number. In dimensional units, this 
corresponds to an initial condition in the region around kf 1.5 x 10 2 mm . Moreover, 
we apply a Gaussian filter in physical space to achieve a beam profile comparable to that 
of the experiment. 

Applicability of the |WT| approach is verified by the calculation of the nonlinear pa- 
rameter J L , equation ( |10[ ), for the numerical simulation, which agrees with the exper- 
iment and is of the order J L ~ 100. Experimentally, we measure the light intensity 
I(x, z) — and not the phases of ip. Therefore the wave action spectrum is not di- 
rectly accessible. Instead, we measure the intensity spectrum, N(k,z) — |/k(z)| 2 , where 
the /c-scaling for in the inverse cascade state is easily obtained from the |KZ| solu- 
tion (75b I and the random phase condition, (see Appendix Q]). This procedure gives an 



intensity spectrum of 

iV k cx£r 1/5 . (121) 
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The wave action spectrum from the numerical simulation is shown in figure [TT] at two 
different distances, we see at z = ram the peak from the initial condition at high wave 
numbers, then at z = 63 mm we see evidence of an inverse cascade, as the majority of 
the wave action is situated towards low wave numbers. However, at low wave numbers 
we do not see the spectrum matching our theoretical KZ prediction (75b I. We showed in 
Subsection 5.13 that the wave action flux, Qk 7 corresponding to the KZ solution (75b I, 
has the incorrect sign for an inverse cascade. Therefore, we noted that the inverse cascade 
spectrum would correspond to a mixed thermal solution with a finite flux contribution. 
This is the probable cause for the lack of agreement with the KZ solution. We observe 
from figure that oc k° at low wave numbers, which may account for the mixed 



thermal solution (as this scaling corresponds to the equipartition of wave action (71b)) 
or for the presence of solitons (c.f. figure [9]). 
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Figure 11: Numerical spectrum of the wave action, at distances 2 = mm and z = 63 mm. 



The experimental and numerical intensity spectra are shown in figures. 12 and 13 
respectively. In both plots, we observe an inverse cascade, resulting in a good agreement 



with the WT prediction for the intensity spectrum (121 1. This agreement may be coin 



cidental, as we have shown that the flux directions are insufficient for a |KZ| wave action 
cascade. However, the intensity spectrum is a non-local quantity, and the system may, 



somehow, conspire to have better scaling properties of the KZ type for the intensity than 
for the wave action. We can only confidently say that we still do not fully understand 
the agreement. 

Verification of the long- wave limit kl^ <C 1 and deviation from integrability is checked 
by considering the ratio of the two nonlinear terms in the LWE| ([7| , which we denote 



12 



and 



as R, and is estimated in Fourier space as R oc k Ij. We observe from figures 
13 that the inverse cascade is approximately in the region k ~ 10 4 — 10 5 m _1 , giving 
an estimation of R ~ 10~ 2 — 1. Note, that if R is too small, then we are close to an 
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integrable system, which would be dominated by solitons and lack cascade dynamics, 
and if R is too high then the |LWE| is a poor approximation for |1D|0WT) 




10 1 10 2 
k [mm -1 ] 

Figure 12: Experimental spectrum of the light intensity, = |/k| 2 at two different distances z. 




Figure 13: Numerical spectrum of the light intensity, = |/kl a t two different distances z. Averaging 
is done over a small finite time window and over ten realizations. 
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The inverse cascade can be observed directly in the physical experiment by inspecting 
the light pattern in the (x, z) plane of the LC cell. Two magnified images of the intensity 
distribution I(x, z) showing the beam evolution during propagation in the experiment are 



displayed in figure 14 For comparison, in figure 14 1 and b, we show the beam evolution 



in the linear and in the weakly nonlinear regimes, respectively. In figure |14fc ,, we set a 
periodic initial condition with a uniform phase and apply no voltage to the |LC| cell. We 
see that the linear propagation is characterized by the periodic recurrence of the pattern 
with the same period, a phase slip occurring at every Talbot distance. This is defined 
by p 2 /A, with p, the period of the initial condition and A, the laser wave length |127j . In 
figure [Tip, we apply a voltage, V — 2.5 V to the LC cell. The initial condition is periodic 
with the same period as in figure |14fc i, but now with random phases. We observe that 
the initial period of the pattern is becoming larger as the light beam propagates along z. 




Figure 14: Intensity distribution I(x, z) showing the beam evolution during propagation; a) linear case, 
b) weakly nonlinear case. 



While the linear propagation in figure 14 1, forms Talbot intensity carpets |128j , with 
the initial intensity distribution reappearing periodically along the propagation direction 
z, the weak nonlinearity in figure |14b , leads to wave interactions, with different spatial 
frequencies mixing and the periodic occurrence of the Talbot carpet being broken. 

In figure 15 we show two intensity profiles taken in the nonlinear case at two different 
stages of the beam propagation. The inverse cascade is accompanied by a smoothing of 
the intensity profile and the amplification of low wave numbers components. 

The intensity distribution I(x, z), showing the beam evolution during propagation 
in the experiment and in the numerical simulation, are displayed in figures [16] and |17| 
respectively. In the high resolution inset of figure [16} we can observe that the typical wave 
length of the waves increases along the beam, which corresponds to an inverse cascade 
process. Furthermore, one can see the formation of coherent solitons out of the random 
initial wave field, such that in the experiment, one strong soliton is dominant at the 
largest distance z. Both the experimental and the numerical results in figures [16] and [17] 
indicate that the total number of solitons reduces. The observed increase of the scale and 
formation of coherent structures represents the condensation of light. Experimentally, 
the condensation into one dominant soliton is well revealed by the intensity profiles I(x) 
taken at different propagation distances, as shown in figure [i~8] for z — 0.3, 4.5 and 7.5 
mm. Note that the amplitude of the final dominant soliton is three orders of magnitude 
larger than the amplitude of the initial periodic modulation. 

Figure 19 displays three experimentally obtained profiles of the PDF of the wave 
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Figure 15: Two intensity profiles I(x) recorded at z = and z = 1.9 mm in the weakly nonlinear regime, 
with V = 2.5 V, showing the smoothing associated with the inverse cascade. 




Figure 16: Experimental results for intensity distribution I(x,z). The area marked by the dashed line 
is shown at a higher resolution (using a larger magnification objectif). 



intensity along the cell at distances z = mm, z = 3 mm and z — 8 mm in lin-log 
coordinates. In a pure Gaussian wave field, we would observe the Rayleigh distribution 
(51) which would correspond to a straight profile of the PDF in lin-log coordinates. 
However, in figure |19| non-Gaussianity is observed with the deviation from the straight 
lines, indicating a slower that exponential decay of the PDF tails. Non-Gaussianity 
corresponds to intermittency of WT Intermittency implies that there is a significantly 



higher occurrence of high intensity structures compared to that predicted by a Gaussian 
wave field and that the system is in the presence of an |MI| process leading to solitons. 
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Figure 17: Numerical results for intensity distribution I(x,z). The frame on the left is a magnified 
section of the initial propagation of the beam. 
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Figure 18: Linear intensity profiles I(x) taken at different propagation distances, z = 0.3, 4.5 and 7.5 
mm. 



We see a similar scenario in figure [20j showing the intensity PDF 5 in the numerical 
simulation. We plot three |PDFfc at distances z — mm, z — 31 mm and z = 63 
mm. The observation of intermittency could be a sign of the development of coherent 
structures (solitons) in the system. 

To allow us to capture the inverse cascade dynamics in more detail and to enable us 
to capture solitons and analyze their interactions with incoherent waves, we perform an 
additional decaying simulation with a lower initial intensity. However, to compensate for 
the lower intensity initial condition, we must run the simulation for much larger times. 
Separation of the random wave and the coherent soliton components can be achieved 
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Figure 19: |PDFfc of the wave intensity within the experimental cell at three different distances along 
the cell, 2 = mm, z = 3 mm and z = 8 mm. Straight lines correspond to the Rayleig h |PDF| s 
corresponding Gaussian wave fields (these fits have the same mean as the respective numerical IPDFE). 




by performing an additional Fourier transform with respect to z over a finite z-window 
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Figure 21: The (k, cj)-plot of the wave field at z = 10 m. The Bogoliubov dispersion relation is shown 
by the dashed green line. 



which results in the 2D Fourier (k, w)-plots of wave action density [45j I49j . 

The numerically obtained (k, w)-plot enables the direct observation of the dispersion 
relation of random waves, and is shown in figure |21| Here, the incoherent wave com- 
ponent is distributed around the wave dispersion relation, which is Bogoliubov-modificd 
by the condensate, (G.4), and is shown by the solid line in figure 21 The distribu- 
tion of the dispersion relation is centered around the theoretical prediction (G.4). We 
observe that the width of the dispersion relation is narrow for large k, a sign of weak 
nonlinearity, and progressively gets wider as it approaches smaller wave numbers. This 
broadening, is an indication that the nonlinearity of the system increases towards smaller 
k, as was theoretically predicted in Subsection 5.12 At wave numbers around zero, we 



see that the theoretical Bogoliubov curve vanishes, corresponding to the region defined 
by equation (105), where the frequency becomes complex. For such wave numbers k, 
|MI| of the wave packets occurs and the |WT| description breaks down. Below the region 
where the Bogoliubov relation becomes complex, we observe slanted lines. Each of these 
lines corresponds to a coherent soliton, whose speed is equal to the slope of the line. To 
examine the time evolution of the system, we have put together successive (k, w)-plots 
into a movie. We observe that the formation of solitons is seen in the movie as straight 
lines peeling with a gradient tangential to the dispersion curve. Moreover, we observe 
the gradual migration of these lines to higher negative frequencies, as the solitons begin 
to grow in size by the absorption of energy from surrounding waves or by merging with 
other solitons. 

Further analysis can be done by separating the wave and the soliton components of 
the (k,w)-plot. We define the soliton region of the (k, w)-plot as the region defined with 
-w m ax/2 < w < 0, (where w k ranges from [-w max , w max ) and w max = *4 ax in the (k,w)- 
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plot). We chose this region to represent the soliton component because we found that 
some wave component of the dispersion relation appears in the extreme negative w-region 
due to the periodicity of the Fourier transform and possibly from an insufficient time 
window for the Fourier transform in z. Consequently, the wave component is defined 
as the region of the (k, w)-plot outside the soliton region. By inverting the Fourier 
transforms of each of the soliton and wave regions separately, we can recover the wave 
field ip(x, z) for each component. This enables us to compare the soliton and wave 
components of the intensity profile, figure [22j and the wave action spectrum, figure [23] 
In figure 22 we plot the intensity profile of light in real space for the soliton and wave 
components. We see that our procedure does indeed isolate the coherent solitons from 
the random wave background and we observe that the main soliton is a least an order of 
magnitude greater than the random wave field. 
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Figure 22: Plot of the intensity profile at z = 10 m of both the soliton and wave components 



In figure 23 we observe that the soliton component of the wave action spectrum, 
is situated towards the low wave number region and scales as oc k . This profile is 
qualitatively similar to the wave action spectrum of the exact soliton of the |1D||NLSE| 
seen in figure [5J The wave component of the wave action spectrum is more widely 
distributed in k-space. Taking into account that the initial spectrum was concentrated 
near k = 1.5 X 10 2 mm~ 1 , we see that, roughly, the inverse cascade range is soliton 
dominated and the direct cascade range is random-wave dominated in this run. We 
also put a slope corresponding to the inverse cascade |KZ| spectrum to show that it 
strongly disagrees with the spectrum in the soliton dominated range. In the direct 
cascade range, the spectrum is consistent with the thermalized spectrum for energy 



(71a |. For comparison, we also plot the direct cascade KZ scaling in figure 23 We see 



that the thermalized spectrum is in a much better agreement with the data than KZ 



This is not surprising, as the decaying simulation is void of any dissipation mechanism, 
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so energy would tend to thermalize. We would expect to see a better test of the |KZ| 
predictions in the forced and dissipated simulations where a non-equilibrium stationary 
state can be achieved. 

Note that in other wave systems, the separation of the coherent and wave components, 
may be a useful technique in observing the KZ scaling otherwise masked by the presence 
of coherent structures. 




Figure 23: The wave action spectrum of both the soliton and wave components of the spectrum, with 
thc|KZ|and thermalized energy predictions for comparison. 



From the same simulation, we plot the [x, z)-plot for the intensity distribution, I(x, z), 
in figure [24) With the longer propagation distance, we observe the formation and evolu- 
tion of solitons from the weakly nonlinear background. In the initial stages of the sim- 
ulation, we observe enlargement of the characteristic wavelength of the random waves 
followed by the formation of numerous solitons out of the random wave background. As 
the simulation progresses, the solitons become more pronounced and begin to behave 
independently from the waves, with an almost random movement through the numerical 
box. A large number of solitons are produced at the beginning of the simulation, but over 
long evolution times, the number of solitons reduces, with the remaining solitons being 
of increasing amplitude. Indeed, there are many soliton merging events. In figure [25b we 
see a magnified section of a merging event from figure [24] At late times, figure [24|shows 
a single dominant soliton in the system that has subsequently grown by absorption of 
energy from fellow solitons and the background wave field. Such a final single-soliton 



state is also observed in the experiment (figure 16 1, albeit at a much shorter propagation 
distance, with the multi-soliton range missing. In the future, it would be interesting 
to set-up a larger experiment which could capture the multi-soliton transient state. In 
figure [25j we plot three zoomed sections of figure [16] to highlight aspects of the soliton 
behavior. In figure [25k, we observe several solitons passing through one another with 
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little deviation in their trajectories. Moreover, the weaker solitons (light gray in color), 
have almost straight trajectories, like freely propagating linear waves, whilst the stronger 
solitons' movement is more erratic. 

In figure |25b , we magnify a merging event between two solitons. The larger soliton 
engulfs the smaller soliton almost without any deviation in its trajectory. Therefore, in 
the numerical simulation, large solitons absorb energy from smaller solitons, reducing the 
total number within the system. In addition, we also observe at z = 11.5 m, a soliton 
bouncing off the larger soliton. As the weaker soliton approaches the larger soliton, it 
slows, before moving away at a fast speed. When two solitons merge and when they 
repel is a key question that still remains to be answered. Finally, in figure |25b , we 
observe two weak solitons propagating together, until around z — 13 to, when they both 
repel each other and disperse back into the random wave field. This shows the break 
up of the coherent structures and the subsequent re-injection of wave action into high 
frequency waves. All these events are parts of the |WTLC| i.e. the cycle of coexistence 
and interactions of the random waves and the solitons and their mutual transformations 
into one another. 

In figures [26] and |27[ we plot the maximum of the intensity and the energy of the 
system with z. To begin, we note that the maximum of the intensity of the wave field is 
always growing, and moreover, we observe that there is a sharp increase in the maximum 
at around z = 11 to. This jump corresponds to the merging event seen in figure [25p . 
Thus, the dominant soliton instantly grows in size once it absorbs the other soliton. This 
merging event is also noticed in figure |27[ where we observe a similar sharp increase to 
the linear and nonlinear energies at z = 11 to, showing that as the dominant soliton's 
amplitude increases, there is a significant increase in the nonlinearity of the system. It 
is interesting that the amplitude of the final dominant soliton continues to grow at the 
late stages of this simulation, e.g. in the range of z between 15 m and 20 m. This 
soliton acts as a kind of "optical vacuum cleaner" which, as it moves, sucks in energy of 
the random waves and remaining soliton debris (i.e. soliton-like structures which are so 
weak that they cannot maintain their coherence for long). 

From the beginning, the nonlinear energy, associated mostly with the first nonlinear 
term in the evolution equation |7]), grows almost by an order of magnitude in size, 
compared to just under double in size for the linear contribution. This indicates that the 
inverse cascade and the subsequent soliton development is associated with an increase 
in the nonlinearity of the system, and hence the breakdown of the |WT| description. 
Moreover, we note that the total energy remains conserved, verifying that our simulation 
is well-resolved. 

In figures 28 we see the soliton and wave component energies evolving with 'time' 
z. The data was made by using the same soliton- wave separation technique as the one 
we described previously, now using a series of (k, w)-plots and computing the energy 
values at each time. In figure [28^ , we observe that over time the total energy of the 
wave component increases, whilst the energy of the soliton component decreases, crosses 
zero and becomes negative. However, this is not indicating a reduction in the soliton 
strength but the opposite, noting that the first nonlinear energy term in the long-wave 
Hamiltonian is negative, this indicates that the soliton component is becoming more 
nonlinear (an increase in the negative energy contribution). This fact is supported by 
figures |28b and [28b, where the linear and nonlinear energies of the soliton component 



are growing in magnitude and becoming of similar values. In figures 28 : and 28 i, we 
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Figure 24: Numerically obtained I(x, z)-plot for a long time simulation with a low intensity initial 
condition to see the inverse cascade and soliton merging. 





Figure 25: We present three close-up images of the (x,z)-p\ot of figure [24] highlighting three types of 
soliton behavior: a) passing through each other, b) merging and c) dissipating. 



have labeled the point in time where two large solitons merge (see figure |25p) . At this 
point, we observe a sharp acceleration in the growth of both the nonlinear energies, which 
corroborates figures 26 and 27 showing a sharp increase in the maximum intensity, linear 
and nonlinear energies in the whole system. Thus figures [28] point towards a system that 
is becoming more nonlinear in time, with energy flowing from waves into solitons. In 
figure [28p we see almost no growth with respect to the linear energy of the solitons. This 
coincides with the analysis of Subsection |5.14[ where the consideration of the conserved 
quantities led us to expect the generation of waves rather than the increase momentum 
of the solitons. This is indeed the figure [28p shows. This is further supported 
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Figure 26: Plot of the maximum of the wave intensity in physical space versus z. 
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Figure 27: Plot of the energies for the long time simulation. We see the conservation of the total 
energy, the growth of the linear and both nonlinear energies (corresponding to the two nonlinear terms 
in equation ||13a|). 



if we inspect figure [24] If we look closely at the paths of the solitons in figure [24] we 
observe that the solitons at late times are not moving 'faster' than those from earlier 
times. Moreover, what is interesting is that in figure |26[ we observe an increase to the 
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maximum intensity of the system when there is only a single soliton present (z > 11 m), 
i.e. the single soliton is still growing in size, but now by the absorption of wave action 
from the linear wave background and not by merging with other solitons. This can be 
viewed as a vacuum cleaner effect, with the soliton 'hoovering' up wave action from waves. 
However, we note that this process does not remove waves from the system, as clearly in 
figure [28) j we observe the wave component of the linear energy is still increasing. Note 
that the growth of the wave energy is consistent with the loss of their wave action when 
the mean wavelength is decreasing. It seems that the 'hoovering' effect is more efficient 
when the soliton is large in size, as opposed to earlier times when solitons are small in 
amplitude. The evidence for this scenario is observed in figures 26 and[28j where, at early 
times, sharp soliton growth (and subsequently the increase of the linear and nonlinear 
energies) is more determined and dominated by discrete mergers between solitons (as 
seen by the sudden spike in growth during the large merger event). However, at late 
times, when no more merging events occur, we still observe a growth in the maximum 
intensity, linear and nonlinear energies but in a more smooth and continuous fashion. At 
asymptotically large times, we expect that the equilibrium state will consist of a single 
stationary soliton, containing the majority of the wave action in the system, and a sea of 
small amplitude short random waves, containing the majority of the linear energy of the 
system. This scenario was originally predicted in [75J H29| and shown using statistical 
mechanics arguments in [551155] . 

In figure 29 we plot the intensity PDF in k-space for two wave numbers, k = 9.8 
mm~~ 1 1 and k = 2.9 x 10 2 mm -1 at z ~ 20 m. The two wave numbers are situated in the 
low wave number (soliton dominated) and high wave number (weak wave dominated) 
regions respectively. We observe at both large and small values of Jk both |PDF| j's 
show signs of depletion of their probability with respect to the Rayleigh distribution. 
However, there is far more depletion observed in the soliton dominated |PDF at k = 9.8 
mm -1 . Probability depletion in the wave dominated scale, k = 2.9 x 10 2 , at the large 
(small) Jk corresponds to a positive (negative) probability flux J 7 ^. We believe that the 
depletion at both large and small Jk's is because of the hoovering effect of the soliton - 
as solitons grow, their hoovering effect (absorption of wave action from waves) becomes 
more apparent and thus more depletion is seen. As the soliton is composed of a wide 
range of wave modes, this effect will occur at a large range of wave numbers (even high 
wave numbers). However, as the solitons merge, grow, or dissipate, they excite high 
frequency waves. This behavior would encourage the transportation of wave action from 
the soliton to small JkS. Indeed, this process would generally only occur at high wave 
numbers and technically would produce an enhancement (negative probability flux) of 
the PDF| tail in figure 29 The balance between the hoovering effect and excitation of 
waves is probably what is occurring at high wave numbers, giving rise to an almost 



(slightly depleted) Rayleigh PDF at k = 2.9 x 10 2 mm -1 . Over time, as solitons merge 
and dissipate, until only one remain, we would observe a reduction in the re-injection of 
wave action at high wave numbers because the mechanisms that allow this (merging and 



PDF 



dissipation of solitons) would cease. At low wave numbers, at k = 9.8 mm , the 
gradually becomes more and more depleted (deviated from Rayleigh) as time progresses. 
As this scale is soliton dominated, we cannot apply our theory of the |PDF| depletion 
in terms of the probability flux JF^ which was developed for weak waves. In fact, such 
a |PDF| corresponds to a certain distribution of solitons in the soliton parameter space 
(A,B) for which we have not got a theory yet. 
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Figure 28: Plot of the evolution of the (a) total energy, (b) linear energy, (c) first nonlinear energy term, 
and (d) second nonlinear energy term, of the soliton and wave components versus z. 



The observed directions of the probability flux directions Jk - toward large Jk's at low 
fc's and to low Jk's at high fc's - can be combined with the picture of the inverse-cascade 
energy flux in the direction of low fc's by considering a 2D flux pattern on the (fc, Jk)- 
plane. This yields a diagram for the |WTLC| presented in figure KM Here, the initial 
condition is at high fc's and low Jk's (weak small-scale random waves). The inverse 
cascade transports wave action, via six-wave mixing, to small fc's and to higher Jk's, 
as the inverse cascade is associated with an increase in nonlinearity of the system. As 
wave action accumulates at large scales, [Ml| kicks in and solitons form, with wave action 
spreading along the soliton spectrum 7ik oc fc°. Note that the fc-modes within the soliton 
are coherent, i.e. correlated with each other. The non-integrability of the system allows 
for solitons to interact with surrounding structures, that enable energy exchange. In 
particular, solitons will collide, occasionally merge or deteriorate via random interactions 
- each of these processes allowing for the emission of energy and wave action to incoherent 
waves, resulting in a reversal of the probability flux T at high fc's. The wave action re- 
injected back into the random wave component fuels the continual process of the inverse 
cascade. This continuous transport of wave action can be seen as a vortex-like flux 
structure in the (fc, Jk)-plane which we call the WTLC in our |1D"1|0WT| system [ID EI]. 
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Note that from figure [29] we see a positive probability flux even at high wave numbers! 
Thus, the WTLC| picture with flux lines not emerging and not terminating anywhere, as 
in figure |30j is only valid for a steady state, whereas figure [29] corresponds to an evolving 
transient where the net flux (produced by the hoovering effect) is still directed from 
waves to (growing) solitons. For the LWE the asymptotical final state is a single soliton 
containing most of wave action and surrounded by almost linear waves containing finite 
energy but almost no wave action. 




Figure 29: |PDF| of intensity J^, V = (<5(sk — ^k)}i at two wave numbers, k = 9.8 mm -1 , and k = 2.9 X 10 2 
mm ■ The black dashed line is the Rayleigh distribution that corresponds to a Gaussian wave field. 
The inset shows the power-law behavior of the PDF 



tails. The black dashed line in the inset is s. 



7.1.2. Forced and Dissipated Simulations 

In this Subsection, we explore numerical simulations of the |LWE| and |SWE| with 
forcing and dissipation. This setup is intended to produce more idealized conditions for 
non-equilibrium stationary states assumed for the realizability of the pure |KZ] solutions, 



( 75 ). Thus, we can test to see if the KZ spectra are realizable in principle, or if they could 
not be form at all, even in the most idealized of settings. In an infinite sized system, both 
the direct- and the inverse-cascade ranges can be realized in the same simulation, with 
forcing acting at an intermediate scale and two inertial ranges either side (see figure [6]). 
However, computational restrictions make it impractical to perform such a simulation. 
Therefore, we simulate each cascade separately, by performing a simulation for each with 
forcing at cither ends of k-space, allowing for just one inertial range. We numerically 
solve the non-dimcnsionalized equation (see Appendix H) with a — 1024. We run the 
simulations with a time step of At = 1 x 10~ 4 in a box of length L = 32ir, with spatial 
resolution of N = 2 12 . The simulations are run so that a non-equilibrium steady state 
is achieved. This is checked by observing stationarity of the total energy % and wave 
action J\f in the system. 
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Figure 30: Diagram of the |WTLC| The high wave number initial condition subsequently gets transferred 
to low k and high Jy by the inverse cascade. At low k, the weak [WT| fails (near the T-junction of the 
dashed lines) and |MI| occurs, solitons form resulting in the wave action getting spread along a <x 
spectrum (horizontal dashed line). The emission of waves from merging or deteriorating solitons re- 
injects wave action at low intensities. 



For the direct cascade setup for the LWE we introduce dissipation profiles at the two 
limits of k-space that removing energy and wave action that is being constantly injected 
by the for cing. We set the low wave number dissipation profile, D^, to be friction, defined 
by (118b I, with the friction coefficient of ^friction = 5x lO^ 1 and apply a hyper-viscosity 
dissipation scheme at high wave numbers with a coefficient of t'hyper = 1 x 10~ 13 . Forcing 
is situated at large scales and is defined by ( 116 ), with amplitude A = 4.1 x 10 4 . Once the 
simulation has reached a statistically non-equilibrium statistical steady state, we analyze 
the statistics by performing time averages over this steady state regime. We have plotted 
the wave action spectrum rik for the direct cascade regime in figure [31] 

We observe that the wave action spectrum in figure [31] does not agree with the 
theoretical KZ prediction (89). There is a clear flat spectrum of the form oc k°. 



This could correspond to several phenomena: solitons, critical balance and a thermalized 
state with equipartition of wave action. However, in figure 32 we have plotted the 
(k, w)-plot for the same simulation at late times. We observe an absence of any soliton 
structures, like the ones we saw in the decaying case (see figure 21). Therefore, we can 
safely reason that the flat spectrum is not associated with the presence of solitons within 
the system. The explanation for the k slope is more likely down to the thermalization 
of the system with respect to wave action around the low to the intermediate wave 
number range. This behavior should be expected due to the lack of agreeing flux signs 
with Fj0rtoft's argument for the |LWE] where a system in such a situation would tend 
towards a modified thermalization spectrum [35] also known as a "warm cascade" [t7] . 
The lack of agreement with |KZ| can be explained by the incorrect flux directions which 
can be most easily understood and determined using the DAM in Subsection |5.13| It 
remains to be understood why we see that the spectrum corresponds more closely with 
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Figure 31: The wave action spectrum nt in a statistically non-equilibrium stationary state for the direct 
cascade simulation of the |LWEl We plot two bla ck lin es to show the theoretical comparisons. The black 
dashed line represent the |WT[ prediction of the |KZ| solution % <x k~ 1 ln(fc£). The black dotted line 
represents the thermalization of wave action. 



the thermalized KE solution with equipartition of wave action, not the one corresponding 
to the equipartition of energy. On the other hand, we observe that the (k, w) -distribution 
is not too narrow around the dispersion curve uj — uj k = k 2 , with the width ~ uik for 
k < 10. Thus, we cannot rule out that the critical balance mechanism could be the true 
reason for the flat spectrum. 

In figures 33 we plot PDF 3 of (a) the wave intensity in the physical space and (b) of 
the Fourier-mode intensity at two different wave numbers. In figure |33[ a), we observe 
an almost perfect agreement of the the [PDF] to the Rayleigh distribution. This indicates 
that the system corresponds to a Gaussian wave field - a good indication that the system 
is in a thermalized state of random waves. Additionally in figure 33 [h) we observe that 
both PDF i are Rayleigh distributed (here the Rayleigh distribution is fitted to the tails 
of the PDF i) . Again, we observe a clear agreement of the tails to the straight line of the 
Rayleigh distribution at large amplitudes. Towards the core of the PDF s in figure 33 'b) 
we see some depletion of the |PDFfr to the Rayleigh distribution the origin of which 
remains unclear. It is possible that this is a signature of appearance of an embryonic 
condensate (the true condensate, solitons, are are not present in this system due to low 
wave number dissipation) - a similar effect was observed in recent simulations of the 3D 



NLSE [47 . 



The inverse cascade displays more unusual behavior. The inverse cascade simulation 
for the |LWE was setup with forcing situated at high wave numbers with a profile given by 
equation (1171 and with low wave number dissipation in terms of friction and high wave 
number dissipation in terms of hyper-viscosity. The parameters used for the simulation 
are A = 7.7 x 10 3 for the forcing amplitude, ^friction = 1 x 10 _1 and ^hypcr = 1 x 10 -11 
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Figure 32: (k, w)-plot for the direct cascade simulation of the |LWE| The green line represents the linear 
wave frequency ui^ = k 2 . 
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Figure 33: (a) Normalized |PDF| of the wave intensity in the direct cascade simulation of the |LWE) The 
straight line corresponds to the Rayleigh distribution, (b) Normalized |PDF| of at two wave numbers 
for the direct cascade in the |LWE] The black dashed line represents the Rayleigh distribution. 



for the dissipation coefficients. As usual, we wait for a non-equilibrium stationary state 
to be achieved before performing any analysis. 

In figure [34] we plot the wave action spectrum from the numerical simulation, time 
averaged once the simulation has reached a non-equilibrium stationary state. For com- 
parison, we also show the KZ prediction for the inverse cascade of wave action. We 
observe a lack of agreement to any power-law behavior to the spectrum. Indeed, we 
see a peak of the spectrum around the forcing region, while the spectrum gradually fall 
off at small and large k. The absence of a |KZ| spectrum can be explained by the in- 
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Figure 34: The wave action spectrum, in a statistically non-equilibrium stationary state for the 
inverse cascade simulation of the |LWE| The straight line represents the |WT| prediction of the |KZ| 
solution of rik <x fc -3 / 5 . 



correct flux directions, but what is really surprising is the omission of any thermalized 
or |CB| spectrum. The |CB| in this case, is identical to the spectrum of the spectrum of 



the thermalization of wave action: oc k . CB has been observed in other optical 



systems, i.e. one described by the 3D NLSE [57]. We found that the spectrum to the 
left of the forcing region seemed to be determined by the level of friction used at low 
wave numbers. By tuning the coefficient of friction, we adjusted the asymptotic value of 
the spectrum at k — 0, and the spectrum would simply gradually decline to that level 
from the forcing region. Dampening of these low wave number modes could prevent 
some non-local behavior that is driving the inverse cascade. Clearly, the dampening of 
these modes is restricting |MI| and the formation of solitons, which may have an essential 
role in the dynamics of the inverse cascade - as we do not observe this behavior in the 
decaying simulation of the previous Subsection. To be certain that we are not observing 
any effects that might be attributed to the presence of solitons, we show the (k, w)-plot 
for the simulation in figure |35| In figure [35] we observe that no solitons are present by 
the lack of any coherent straight lines around dispersion curve. However, if one looks 
closely, we do observe the partial formation of a soliton situated on the dispersion curve 
at k = 15. However, with the presence of dissipation this soliton does not 'peel' from 
the wave component. Moreover, inspection of the intensity distribution (not shown), 
does not provide any indication of coherent soliton like structures distinguishable from 
random waves. The (k, w)-plot is also in good agreement to the linear dispersion relation 
u>k = k 2 shown by the green dashed curve, but the width of the plot is as wide as uik in 
the inverse cascade range. This implies that the system is indeed in a strongly nonlinear 
regime which could point at CB but, as we have already mentioned, the |CB| spectrum is 
not observed. We have also done simulations with weaker forcing amplitudes for which 
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|WT| was weak, and the result was similar to the one reported here. In fact, the reason 
that we increased the forcing was an attempt to strengthen the wave interactions in order 
to "push" the inverse cascade through to the large scales. 



1500 



1000 



S 500 







-500 




-40 -30 -20 -10 




k 



10 20 30 40 



Figure 35: (k, aj)-plot for the inverse cascade simulation of the |LWE| The green line represents the linear 
wave frequency oj^ = fc 2 . 



In figures 36 we plot the PDF, of (a) the wave intensity in the physical space and 



(b) of the Fourier-mode intensity at two different wave numbers. Similarly to the direct 
cascade, we observe agreement to the Rayleigh distribution of a Gaussian wave field. The 



PDF of the wave intensity in figure 36 'a) is in perfect agreement to the unfitted Rayleigh 
distribution, whilst we observe good Rayleigh agreement of the tails of the |PDFfc in 
figure 36 'b) at two wave number either end of the inertial range. Moreover, we observe 
a depletion of the probability at low wave mode intensities at both k = 2 and k = 10, 
the same phenomenon as seen in the direct cascade. 

In the simulation of the forced/dissipated inverse cascade, we see no agreement of 
the |KZ| thermal or |CB| spectra. Moreover, we observe that the spectrum seems to be 
determined by the value of the lowest mode. This suggests that the inverse cascade 
may be associated to some non-local effect established by |MI| Indeed, we find the 
forced/dissipated simulation to provide more questions than answers to the behavior of 
|lD||OWTl and clearly more work in this area is needed to fully explain these results. 

7. 2. The Short- Wave Equation 



We investigate the SWE ^ by numerical simulations of the non-dimensionalized 
model of Appendix | H| We apply the same forcing and dissipation profiles that are 
described in Section [6j We dissipate at low wave numbers using a hypo- viscosity profile, 
while using hyper-viscosity at high wave numbers. 

We begin by discussing the direct cascade simulation. We force the system with the 
profile (116), with amplitude A = 1.6 X 10 2 and use the dissipation rates of Vhypo = 
10 , t'hypor = 1 x 1CT 11 for the direct simulation. 

61 



1 x 




Figure 36: (a) Normalized |PDF| of the wave intensity in the inverse cascade simulation of the |LWE| The 
straight line corresponds to the Raylcigh distribution, (b) Normalized |PDF| of Jk at two wave numbers 
for the inverse cascade in the |LWE] The black dashed line represents the fitted Rayleigh distribution. 



In figure 37 we plot the wave action spectrum for the direct cascade averaged over 
a time window once the system has reached a statistically non-equilibrium steady state. 
We observe a good agreement to the WT prediction (76a). However, notice that there is 
some slight deviation of the spectrum in the middle of the inertial range. This is likely to 
be caused by an insufficient time average of the statistics. The agreement is surprising, 
as we showed that the KZ solution, ( 76a I , does not produce convergence of the collision 
integral in the UV limit. However, we may be observing a non-local spectrum that is 
close to the local WT| prediction. 




Figure 37: The wave action spectrum % in a statistically non-equilibrium stationary state for the direct 
cascade simulation of the |SWE| The straight line represents the |WT| prediction of the |KZ| solution of 
rik oc fc 7 / 5 . 



62 




Figure 38: (a) Normalized |PDF| of the wave intensity in the direct cascade simulation of the |SWE| The 
straight line corresponds to the fitted Rayleigh distribution, (b) Normalized |PDF| of Jk at two wave 
numbers for the direct cascade in the |SWE] The black dashed line represents the Rayleigh distribution. 



In figures [38] we plot the normalized [PDF for the wave intensity in the x- and k-spaces 
respectively. We observe in figure [38^ a) that there is a slight deviation to the Rayleigh 
distribution (corresponding to a Gaussian wave field) at high intensities. This implies 
that there is a larger than expected occurrence of high intensity structures - a sign of 
wave intermittency. The |S WE] does not seem to produce solitons, but we can numerically 
check this fact my investigating the (k, w)-plot in figure 39 WT theory implies that such 
behavior should exist as part of a wave breaking process. In addition, in figure 38 a) 
we see a large accumulation of the probability at low intensities, far more than what is 
predicted by the Rayleigh distribution indicating an enhanced presence of dark regions 
void of light. On the other hand, figure |38^b) is the PDF] of Jk at two wave numbers 
at either end of the inertial range. We see that the PDF| at k 
close to the Rayleigh distribution 



1 remains relatively 
with only a minor depletion in the tail. At k = 7 



we observe a large enhancement of the |PDF| at large Jk, corresponding to a negative 
probability flux J-^. Therefore, we see behavior that is similar to that which is described 



by the WTLC 30 Observation of the (k, cj)-plot, figure 39 shows that actually there 
are no coherent structures present in the system. This conclusion is made by the lack 
structures at negative uj. Moreover, we see a clear accumulation of the dispersion curve 
situated around the linear wave frequency uik — k 2 . This indicates that the system is in 
the weak nonlinear limit. 

The inverse cascade simulation of thc |SWE] is produced with the following parameters: 
The forcing amplitude is set to be A — 4.8 x 10 2 , while the dissipation rates are given 
by t'hypo = 1 x 10°, ^hypcr = 1 x 10~ 8 . Similarly to the direct simulation, we wait until a 
non-equilibrium stationary state is achieved, and then perform averages on the statistics. 

In figure |40[ we observe a good agreement, for almost a decade in k-space, with 
the KZ scaling (76b). At low wave numbers, we observe a slight accumulation of wave 
action before the dissipation occurs. However, analysis of the locality, showed us that 
the inverse |KZ| solution is invalid - similar puzzling behavior as in the case of the direct 



cascade simulation of the ISWEI 
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Figure 39: (k, w)-plot for the direct cascade simulation of the |SWE| The dashed green line represents 
the linear wave frequency u>y. = k' 2 . 




Figure 40: The wave action spectrum in a statistically non-equilibrium stationary state for the inverse 
cascade simulation for the |SWE| The straight line represents the |WT| prediction of the |KZ| solution of 
n k <x fc" 9 / 5 . 



In figures |4l] we plot the normalized |PDF| for the wave intensity in the physical 
space and for the Fourier-space intensities Jk- In both plots for comparison we insert 



(>4 




(a) 



Figure 41: (a) Normalized |PDF| of the wave intensity in the inverse cascade simulation of the |SWE| 
The straight line corresponds to a fitted Rayleigh distribution, (b) Normalized |PDF] of Jk at two wave 
numbers for the inverse cascade in the |SWE] The black dashed line represents the Rayleigh distribution. 



lines corresponding to a Gaussian wave field. We observe in figure 41 'a) that there is a 



clear deviation at low intensities from Gaussianity. This gives a clear indication of WT 
intermittency, with the enhanced frequency of occurrence of dark spots - far greater than 
what is predicted by a general Gaussian wave field. For J k seen in figure [41] we see a 
close agreement at intermediate Jk with the Rayleigh distribution given by the black 



dashed line, but at large Jk we see some slight deviation. This is a possible sign of WT 



intermittency, with the depletion of the PDF tail from Gaussianity. Inspection of the 
(k, w)-plots in figure 42 shows a wide spreading of the wave distribution around the linear 
dispersion curve uJk = k 2 indicating a nonlinearity level that is quite strong. Moreover, 
we observe asymmetric behavior of the wave distribution that clearly shifts from negative 
to positive wave numbers in figures 42 b) and 42 d) . This could indicate that the there 
is some symmetry breaking instability associated with the SWE Furthermore, there 



is an interesting 'hair'-like structure around the dispersion curve in figures [42] which 
align together in a certain direction. It consists of lines in the (k, w)-plot indicating a 
set of coherent structures propagating with speeds that are determined by their slopes. 
The 'hair'-like lines are not straight which indicates that the speed of such 'solitons' 
varies in time. Also note that these structures never 'peel off' from the wave dispersion 
curve like it was the case in |LWE and, therefore, they never de-couple from the wave 
component. In figure |43|a) we plot the distribution of intensity at late times of the 
simulation over the whole domain. We observe strange behavior, in which soliton like 
structures are seen bundled together and weakly interacting, while propagating as a 
bundle with constant speed. These structures are clearly the hair-like objects seen in 
the (k, oj)-plots of figures 42 If we take a closer inspection of the coherent structures in 
figure [43] we see that they meander around the center of the bundle, but remain confined 
to the bundle. We don't yet know how to explain this behavior theoretically or how it is 
related to the KZ spectrum. In figure 43 ^b) we show a time slice of the intensity /. We 
see a low wave number coherent structure with high frequency peaks propagating upon 
it. 
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Figure 42: (k, oj)-plots for the inverse cascade simulation of the |SWE| at four different times r = 13.67, 
18.65, 28.60 and 33.57, where r is in units of the linear evolution time Tl of the slowest wave mode in 
the system. The dashed green line represents the linear wave frequency uj^ = k 2 . We observe asymmetry 
of the dispersion curve in figures (b) and (d) with respect to positive/negative wave numbers. 



8. Conclusions 



We introduced the field of OWT where waves of light propagating through a nonlinear 
medium weakly interact. We presented the first ever experiment of [OWTj which is 
realized in a |f D| setup based on a |LC| cell, and that serves as a basic example to which we 
refer when developing the theoretical framework for |lD|[OWT| We showed, by starting 
with an initial condition of weakly nonlinear waves, we were able to, by the inverse wave 
action cascade described by |WT| theory, to show the formation of strongly nonlinear 



solitons. One intriguing aspect of OWT is the inverse cascade, where wave action is 
transported to large scales. This process is analogous to |BEC| in super-cooled alkali gases. 
We demonstrated this process in the experimental implementation. We showed that by 
injecting high frequency waves into a nonlinear |LC| cell, we observed a wave turbulent 
inverse cascade, followed by the condensation of incoherent photons into coherent solitons. 
Moreover, we supported the experimental findings by performing numerical simulation 
of the |LWE| Analysis of the inverse cascade by |WT| theory predicted that a pure |KZ 
state was not realizable. Although the |KZ| spectra are unrealizable, the |OWT| system 
is still a very interesting object which is reach in physical effects and mathematical 
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Figure 43: (a) Distribution of I(x,z) at late times in the simulation of the |SWE) (b) We plot a time 
slice I(x) at distance z = 50001. In (a) we observe coherent structures propagating with constant speed. 
In (b) we observe a snapshot of the condensate structure. 



structure. We show how coherent solitons are formed by nonlinear interactions out of 
random incoherent waves, and analyze the importance of the soliton behavior with the 
random wave background. Additionally, we determined that the inverse cascade can be 
described by a mixed wave action spectrum of thermal and non-equilibrium parts. The 
inverse cascade is associated with the growth of nonlinearity of the system, gradually 
resulting in the formation of solitons by |MI| We observed solitons in both the experiment 
and numerical simulations, and found that the system relaxes to a state of a single 
dominant soliton. During the transition to this final state, we observed interactions 
between solitons and the random wave background. With the aid of WT| theory, we 
developed a strong WT description for the behavior of the system in the form of a 
|WTLC| The |WTLC diagrammatically represents the coexistence of coherent solitons 
and the incoherent random wave background. 

In addition to the experimental setup, we performed numerical simulations of the 



LWE and SWE in non-equilibrium stationary regimes in forced-dissipated systems. This 



allowed for the subsequent testing of the KZ solutions. Theoretically we discovered that 
the |LWE"| emits fluxes of opposite sign to what is necessary for the |KZ| states, and this was 
further verified by a lack of agreement of the wave action spectrum to the |KZ| solutions. 
Conversely, this problem does not occur for the |SWE| However, we showed that both 
the IKZI for the ISWEI are non-local and thus should be non-realizable. Nevertheless, 



we observed good agreements to the KZ predictions for both the direct and inverse 
cascades. Further investigation into the fluxes yielded a non-constant energy flux for the 
direct cascade - a constant flux being an essential requirement for a local |KZ| solution. 
This indicates that we observed a possible non-local wave action spectrum and if so, we 
are led to ask what non-local process determines this power law? The answer to this 
can be determined with a development of a non- local theory for [TD"1|OWT| in th e short- 
wave limit - which is a future goal of ours. The inverse cascade for the [SWE yielded 
an excellent agreement to the |KZ| solution and a good indication of a constant negative 
flux of wave action. This suggests that the inverse cascade is of [KZ] type, contradicting 
the non-locality analysis and thus raising questions in itself. Thus further analysis of the 



G7 



numerical prefactor of the wave action spectrum is required and a possible explanation 
of the observations via a non-local theory. Further, in the inverse cascade simulation of 
the |SWE[ we found a new fascinating physical phenomenon - generation of meandering 
solitons organized in crisscrossing bundles. Theoretical explanation for such a behavior 
or possible experimental implementations are yet to be constructed. 

In trying to tackle |OWT| we have ended with more questions than answers. We 
have observed very reach behavior with numerous interesting physical effects which po- 
tentially could come useful for technological applications insofar as optical solitons are 
considered as potential candidates for optical interconnectors and information process- 
ing. However, there remains a lack of a solid theory for explaining the observed results, 
and some of these results appear to be in contradiction with the classical WT| analysis. 



Thus, we consider this review as a report summarizing the first stage [OWT studies, with 



an open discussion of the observed results, theoretical advances and remaining puzzles. 
Our analysis has raised several new questions that need to be addressed with additional 
experiments and numerical simulations. There is yet to be done an experimental in- 
vestigation of higher dimensional |OWT[ which we believe is a natural extension to the 
research presented in this report. In addition, the description of the evolution of solitons, 
with the possible understanding of their behavior and manipulation would be valuable to 
any future applications in nonlinear optics. In this context [OWT| represents a challenging 
step forward, providing the occurrence of the spontaneous increase of the coherence of 
light through weakly nonlinear wave interactions. 
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Appendices 



A. The Canonical Transformation 

In this Appendix, we will outline the details of applying thc |CT| to a six- wave system 
with non-resonant four- wave interactions. This procedure is adapted from [3], where a 



CT was used for converting a non-resonant three-wave system into one with resonant 



four-wave interactions. The fundamental goal of the CT is to remove the leading order 
non-resonant four-wave interactions by re-defining the canonical variable by a quasi- 
nonlinear canonical transformation. This will preserve the system's linear dynamics 
whilst removing the lowest order non-resonant dynamics. The result will be a system 
that contains resonant wave interactions at its lowest order. 

To begin, we must define a new canonical interaction variable, Ck, such that it sat- 
isfies an auxiliary Hamiltonian, T-L a ux, in the sense that it contains arbitrary interaction 
coefficients at all orders^] We can define a general auxiliary Hamiltonian, up to order 
six, for Ck as: 



tin 



\ E S 1 / (ctc 2 c* 3 + c.c.) + 1 £ 8^ (c lC2 c 3 + c.c.) 

1,2,3 1,2,3 

i 1 if>l,2 rl,2 * * 1 \^ -p-1,2,3 rl,2.3 / * , \ 

f 4 Z^ ^3,4 <>3,4 C 1 C 2C 3 C 4 + - 2^ X 4 <V (CiC 2 C 3 C 4 + C.C.) 

1,2,3,4 1,2,3,4 

^ E y 1 ' 2 ' 3 ^ 1 ' W (ciC 2 C3C 4 + c.c.) 

1,2,3,4 

^ E i 1 ' 2 ' 3A5 '5 1 ^ 4 ' 5 (c lC2 C3C4C 5 +C.C.) 



1,2,3,4,5 



1 61.2,3.4 rl.2,3,4 / * , \ 

+ ^ Z^ B 5 S 5 (C1C 2 C 3 C4C 5 +C.C.) 

1,2,3,4,5 

+ — 2^ C 4,5 $4,5 (CiC 2 C 3 C 4 C 5 + C.C.) 
1,2,3,4,5 

1 tt rl ,2 ,3 rl,2.3 * * * 

+ ^ Z^ ^4,5,6 d 4,5,6 ClC 2 C 3 C 4 C 5 C 6 

1,2,3,4,5,6 

1 ^ Al. 2,3,4,5 rl. 2,3,4,5 / * . \ 

+ 7^ Z^ ^6 ' ' <V (ciC 2 C 3 C 4 C 5 C 6 + C.C.) 

1,2,3,4,5,6 

1 61.2,3.4 cl. 2,3,4 / * * . n 

+ ^ Z^ ^5,6 5 5,6 (C1C 2 C 3 C 4 C 5 C 6 + C.C.) 

1,2,3,4,5,6 

• E S 1AM " 6 ' 6 « 1AM ' 8 'nc lCa c 3C4CB c 8 + c.c.). (A.l) 



1,2,3,4,5,6 



4 We have defined the auxiliary Hamiltonian up to six-wave interactions, which is sufficient for our 
systems. This will be verified a posteriori. 
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The auxiliary Hamiltonian is completely arbitrary and as such we denote all interaction 
coefficients with tildes to emphasize this. The idea is to construct a new Hamiltonian, 



( A. 1 1 , for our system in terms of Ck using the original Hamiltonian ( 16 ) to determine cer 



tain auxiliary interaction coefficients, whilst allowing us to explicitly choose the remaining 



of (161 



arbitrary interaction coefficients in ( A.l I to eliminate the non-resonant interaction terms 



The |CT| is weakly nonlinear, thus preserving the linear dynamics of the system. We 
utilize the fact that the time evolution operator is canonical 3], and then use the Taylor 
expansion of ak around a(k, 0) = c(k, 0) giving 



a(k, z) — c(k, 0) + z 



9c(k, z) 

dz 



z" i <9 c(k, z) 
2 



(A.2) 



auxiliary Hamiltonian, (A.l), i.e 



The coefficients of the CT in (A.2) can be computed using relation (20) applied to the 



d 2 c k 
dz 2 



.SH a ax 
. d 5H m 

! 

dz Set 



(A.3a) 



(A.3b) 



These coefficients will be of the form of summations, involving the canonical variable 
Ck and the interaction coefficients of the auxiliary Hamiltonian. Coefficient (A. 3a I is 



expressed below, and it simply determined by considering the variational derivative of 
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the auxiliary Hamiltonian with respect to c k : 



( dc(k,t) 
\ dz 



1.2 



ck.l * 
2 C X C2 



+ 2 E {U^Y S^clc; + \ £ 8%c lC2 c 

1,2 1,2,3 

, 1 v-l,2,3xl,2,3 , o / ^k.l,2\ * ck,1.2 * 

+ 6 k °k CiC 2 e 3 +3(Xj j S 3 Cl 

1,2,3 

+^ E ( v k : j ; )'^ k : - v, 



c 2 c 3 



1,2,3 



24 



E f^ 1 



,2,3,4 



2 L 3 



rk,l,2,3,4 * * * * 
' L 1 L 2 L 3 L 4 



1,2,3,4 



1 51,2,3.4 rl, 2, 3,4 . , /5k, 1,2, 

+ ^4~ Z^ 5 k Ci €20304 4-4^4 ' 

1,2,3,4 

4-— 2^ 26 k,4 \,4 cic 2 c 3 c 4 + 3 I C 3 ; 4 
1,2,3,4 

, 1 4^1.2,3 s-1,2,3 * * 

+ ^ T k,4,5 <5 k,4,5 C l C 2 c 3C 4 C 5 

1,2,3,4,5 

+4 E « 



ck.l, 2, 3 * * * 
O4 c l c 2 c 3 c i 



ck,l,2 * * 

dg 4 CJC2C3C4 



,2,3,4,5 cl, 2,3,4,5 



ClC2C 3 C 4 C 5 



1,2,3,4,5 
5( Q k ,l,2,3,4 



ck,1.2,3,4 
5 5 



C 1 C 2 C 3 C 4 C 5 



1 r, 51,2,3,4 rl, 2,3,4 * 

■43 Z^ 2R K5 °K5 CXC2C3C4C5 



1,2,3,4,5 



-4 i? 



" ? k.l,2,3 



srk,l,2,3 * * * 
"4,5 c l c 2 c 3 c 4 c 5 



I ^ ^gk, 1,2,3,4,5^ * jk.l. 



2,3,4,5 * * * * * 
C 1 C 2 C 3 C 4 C 5 



1,2,3,4,5 



(A.4) 



The second coefficient, (A.3b|, is derived by taking the time derivative of (A.4) and 
then re-substituting the expression of Ck from (A.4). As one can expect, this leads to a 
lengthy formula which we omit from the review due to space restrictions. With the full 
expression of the CT found in terms of CkS, we then substitute the CT (A. 2) into the 
Hamiltonian (16). The result will be an expression of the Hamiltonian in terms of the 
new interaction representation variable Ck- The Hamiltonian will involve the interaction 
coefficients from the auxiliary Hamiltonian (containing tildes) and T 3 ' 4 ■ The aim is to 
eliminate the non-resonant four-wave interaction contribution by selecting the values of 
the arbitrary interaction coefficients. Moreover, as the system conserves wave action, we 
find that the |CT| does not result in any additional odd A- wave contributions, hence all 
arbitrary interaction coefficients for odd orders are automatically zero. 

71 



We find that elimination of all non-resonant four-wave contributions can be achieved 
by selecting 

-*(l* 9 )* 



T>1,2 _ 
J 3.4 — 



(A.5) 



Relation (A.5) is valid as the denominator does not vanish due to a lack of resonant 



four-wave interactions (i.e. no non-trivial solutions to the four-wave resonance condition 
(25)). After selecting (A.5), we find that Hamiltonian (16) reduces to 



1 

36 



E 

1,2,3,4,5,6 



1,2,3 
4,5,6 



i (u>x + UJ 2 + W 3 — UJ 4 - 0J 5 - Ul e ) 



xW, 



1,2,3 
4,5,6 



(A.6) 



Hamiltonian (A.6) represents Hamiltonian (16) in the new canonical variable Ck, up to 



the leading resonant wave interaction, in this case being of order six. Notice, that within 



rl,2,3 



the six- wave contribution there is still an arbitrary contribution, W 4 ' 5 ' 6 , 
the auxiliary Hamiltonian (A.l). However, its prefactor, lji + uji + LJ3 



that arises from 
- W4 — W5 - uj e . 



vanishes when the six-wave resonance condition 



ki + k 2 + k 3 = k 4 + k 5 + k 6 , 



(A. 7a) 



w(ki) + w(k 2 ) + w(ka) = w(k 4 ) + w(k 5 ) + w(k 6 ), 



(A.7b) 



is satisfied. Ultimately, W^'g g does not provide a contribution to the nonlinear wave 
dynamics. Therefore, we may set 5 ' 6 to equal anything without altering the nonlinear 

""123 • • 123 

wave dynamics. For convenience we select W 4 ' 5 ' 6 to equal minus the difference of W 4 ' 5 ' 6 



from the value taken when the resonance condition, ( A.7), is satisfied. For instance, if we 



decompose W^'^'g into two parts, the first being its value when the resonant condition is 



"4,5,6 



satisfied, say rc wI^'q and the second being its residual 
i.e. we can then express 



yv A,5,6 — KV 4,5,6 



flC W l,2,3 



4,5,6 ) 



1,2,3 
4,5,6 



flCW.2,3 



By choosing 



W 4,5,6 



-L W 1,2,3 



,1,2,3 
4,5,6 



(A. 



(A.9) 



the arbitrary interaction coefficient, W^'g, directly cancels with the residual contribu- 
tion ^Wlfl- This means that the six-wave interaction coefficient will equal its value on 
the six-wave resonance condition over the whole of the domain. The final result yields 



the Hamiltonian expressed as equation ( 26 1 



B. Phase and Amplitude Averaging 



In this Appendix, we will perform phase and amplitude averaging over a RPA field 
for t he derivation of the evolution equation for the|GF| To achieve this, we express the 



GF 



in powers of e, and represent each order in terms of the variable oi . At this stage, 
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we can average over phases by assuming a |RPA| field. Once, phase averaged, we can 
average over amplitudes and compute the evolution equation. 



To begin, we substitute the e-expansion of 6k, (36 1, into definition (43 1 and consider 
terms up to 0(e 8 ). Further simplification can be achieved by Taylor expansion of each 
exponential thus obtaining 



,(A)A k i<ri 2 



(£) A *(^ 1> + ^ (0) 



OL 1 



)|2 



+ 



2 f 



6 (0)* & (D 



« 2) 



+ 0(e 12 ). 



(B.l) 



We now proceed by averaging each term in expansion (B.l I. The leading Oie ) contri- 
bution of Zk is given by 



r(0) 



J (^)^I"L° , I : 



(0) 



2 k (0), 



where J, 



(0) 



|6[°^ I 2 is the amplitude of the interaction variable b^K 



(B.2) 



For clarity, we 



shall omit the superscript of which will thus become JkQ The leading order of the 
e-expansion is labeled as -2-k(0), corresponding to the initial value of the GF The next 
order yields 



L \ 



r(l) 

-k \ \ 27T J ' >k ' " l"k "k "k "k 

and subsequently at the final order to be considered, we obtain 



(B.3) 



r(2) 



L 

L 

2^ 



Ak 

2 



2Jki5k i) i 2 + ( fek o) e) 2 + (e^ 



(B.4) 



Having expressed 2k in terms of 6k, we now proceed by evaluating the ensemble averages. 
We must take the phase average before considering the amplitude average. In the 

1,(1) |2\ 



expansion for the 

h (.0) b (2) 



GF 



we can identify three distinct terms, these are (6 k u k 



l(0).(1)* 



and (b^'b 1 ^ 1 *). We implement Wick's rule, which states that the only contributions are 
those where one can pairwise match the wave numbers. However, due to the trivial 
pairings leading to the nonlinear frequency correction, we are restricted to which wave 
numbers can be matchecJ^J Consequently, for (b^b^*) no pairings can occur and thus 

the average is zero. This is also the case for ((b^b^*^j )<P an< ^ so gi yes zero. 



5 This will not alter the the accuracy of the derivation. 

6 Note that now is not the amplitude of the full interaction variable b^, as defined in Subsection 5.1 
7 See the definition of set /C 4 ' 5 g. 
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In the case of (|& k | 2 ), we are able to match the wave numbers in such a way as to 
not contradict the restrictions of the summations. Applying Wick's rule, we find that 

x [Si (Mi + <Mo) + #0 (<^ii + <M) + tfi (S 5 9 S% + S 5 w 6l)] ■ (B.5) 

This gives 



^ E IW4 k 5 2 6 3 l 2 ^J 3 J4J 5 J6 # B % 3 |A T (u, 4 fc ; 5 2 ; 6 3 )| 2 . (B.6) 



/C 4,5,6 



Applying a similar approach to (& k b k 2 ^* 



we arrive at 



(4 O) ^ 2> )0 = ^ E l W 4,5,6| 2j k ( J 2 + J3) JaJbJsS^Et 

r 2,3 
4,5,6 

-4 E |W4 k 5 2 6 3 | 2 ^kJ 2 J3(J4J5 + J4J6 + J 5 J 6 )^. 



fc.2,3 



k,2,3 
4,5,6 



4,5,6 

,-* / 4.5,6 fe,2,3 

x_ft T lu kA3 ,u 4i5)6 



(B.7) 



By defining the wave action spectrum n k , and the derivative of the |GF| with respect 
to Ak as 

L 
2n 



n k = 



k/J> 



(B.8a) 



5A k 



) J k e(^) AkJk 



(B.8b) 



and using our formulae (B.6) and (B.7), equation (44) can be written in the form, 
Z k (T)-Z k (0) = 



2tt 



4 f! 
12 



E iv* 



2,3i2 rk,2,3 
5,6 I °4,5,6 



A k 2 k + A; 



<9Z k 



9A k 



+2A k 

-2A k 
x3? 



9A k 
9A k 



A ( fe,2,3\ 
A T (w 4 ,5,6 J 

(n 2 + n 3 ) n 4 n 5 n 6 5R 



En 



fe,2,3 



fc,2,3 



4,5,6 J ' W 4,5,6 



™2™3 («4«5 + n 4 n 6 + n 5 7l 6 ) 



TP* / 4,5.6 fc.2,3 
T [ w fc,2,3' W 4,5,6 



(B.9) 



up to C(e 12 ) corrections. 

In equation (B.9), we now take the large box and long evolution time (or equivalently 
the weak nonlinearity) limits. The order in which the limits are taken is essential. We 
must take the large box limit, N — > 00, before the weak nonlinearity limit, T ~ 1/e — > 
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oo, otherwise the width of the nonlinear resonance broadening will be smaller than the 
frequency grid spacing, and will therefore result in a zero contribution. For the large box 
limit, the summations are replaced by integrations of the form: lini/v_>. 00 ^2 3 4 5 6 ~ 

( in ) 5 / ^2 dk.3 cDs.4 dk.5 dk.6 and the Kronecker delta function is transformed into a 
Dirac delta function with a prefactor of 2n/L. 

In the weak nonlinearity limit, we can calculate the asymptotic limits of At{x): 



and of Et(x, y): 



lim |A r (a;)|' 

T— >oo 



lim mE T (-x,x) 

T-s-oo 



2ttTS(x) 



ttTS(x) 



(B.10) 



(B.ll) 



Moreover, as the intermediate time T, at which point we want to seek a solution, is 
smaller than the nonlinear evolution time Tnl, we can make an approximation for the 
time derivative of by: 

Z k (T) - Z k (0) 



2 k 



T 



(B.12) 



The large box and weak nonlinearity limits, together with results (B.10), (B.ll I and 



(B.12 ), gives the continuous description of the evolution of the GF represented in equation 
(45T 



C. The Zakharov Transform 

The|ZT|is used to find steady state solutions to thc|KE| The transformation utilizes 



symmetries of the KE to overlap six disjoint regions of the domain of the collision integral 
on to one another. From the resulting integrand, one can straightforwardly derive the 
steady states of the |KE| In Subsection |5.9[ we presented the transform on one of these 
regions. In this Appendix, we give the four remaining transformations given by 

~ kk-i ~ k ~ kkq ~ kk a ~ kkc , , 

k x = — l -, k 2 = — , k 3 = ^, k 4 = ^ 7 i and k 5 = ^, (C.la) 
k 2 k 2 k 2 k 2 k 2 

with Jacobian J = — ^k/k 2 ^ . 

k 1 = ^i, k 2 = ^, k 3 = ^, k 4 = ^ and k 5 = ^, (C.lb) 
k 3 k 3 k 3 k 3 k 3 

with Jacobian J = — ^k/k 3 ^ . 

ki = ^, k 2 = ^, k 3 = ^, k 4 = ^ and k 5 = 5^-, (C.lc) 

^4 ^4 ^-4 ^4 ^4 



with Jacobian J = — I k/k 



with Jacobian J = — (^k/k 

We re-iterate that these transformations are only applicable if the collision integral 
is proven to be convergent upon the KZ solutions. If not, then the KZ solutions deemed 
non-local and the solutions found using the ZT are unphysical. 



D. Expansion of the Long- Wave Six- Wave Interaction Coefficient 



As the LWE| contains two nonlinear terms, the resulting four- wave interaction coef- 
ficient (18a), will also contain two expressions. WT theory relies on that fact that the 



interaction coefficients are scale invariant. However, both of these expressions express 
different scalings with respect to the wave number k (where all four wave numbers are of 
the same order, i.e. fci, &2, fe, foi cx k). However, for the LWE| this does not pose a prob- 
lem because the long- wave limit kl^ <C 1, provides a small parameter for us to determine 
a leading order, scale invariant, contribution to the four- w ave interaction coefficient. The 

can be represented as a sum of 



final six- wave interaction coefficient VV^'g'g, of the 
three contributions: 



LWE 



1,2,3 
4,5,6 



1 W 1,2,3 
rV 4,5,6 



1,2,3 
4,5,6 



1,2,3 
4,5,6 ' 



(D.la) 



where each term is of a different scaling with k. The exact expressions are given by 
3 6 



1 



E E 

ij,m—l p,q,r=4 



lrpp+q-i,i lrpj+m-r,r 



lrpi+j-p,p lrpq+r 

i,j q,r 



>jj- 



j-\-m—r,r 



q-\-r—m,m 
q,r 



(D.lb) 



4.5,6 



and 



1,2,3 
4,5,6 



3 

E 

i : j,m— 1 4 
i^j^m^i p^q^r^p 

3 6 

E 

i,j,m=l p,g,r=4 
i^j^m^iip^q^r^p 



lrpp+q-i.i 2rpj+m-r,r 



j+m—r,r 
■ 



6 2rpp+q-i,i lrpj+m-r,r 
p,q 



■ 

J, 771 



lrpi+j-p,p 2rpq + J 771,771 

i,j q,r 

q + r — 777,777 

LiJq,r 

2jii+j-p,P lrpq + i 771,771 

" "' (I).lcj 



g + r — 771,777 

q,r 



3 6 2rpp+q~i,i 2jij+rn-r,r 2rpi+j-p,p 2rpq+r-m,m 

y y P , q j,m ^3 q± (Dld) 



i, j, m—1 p,q,r—A 
i^j^m^ip^q^r^p 



UJ 



j-\-m—r,r 



q+r—m,m 
<Jq,r 



where 1 T and 2 T are defined in equation (18a I. 

The three contributions of l wI'^q given in equations (D.l), each are of different 



1,2,3 
4,5,6 



orders in kl% - the small parameter arising from the long-wave limit. Therefore, W, 
mathcalO ((fcZ^) 4 ) is the leading contribution to the six-wave dynamics, followed by 
2 W 



1,2,3 
4,5,6 



O ((A^) 6 ) and finally 3 W, 



1,2,3 
4.5,6 



O ((kl^ 



The first term in expansion (D.l), 1 W4' 5 ' 6 , is generated from the coupling of the 
leading cubic nonlinear term ( 1 T) of |7| with itself. If we were to consider thc |lD|NLSE 
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then the six-wave interaction coefficient would only consist of this contribution ('""WVg'g), 

should imply that the six-wave interaction 



3,4 



0. Integrability of the 



ID 



NLSE 



coefficient 1 W 4 '5'g is zero when the six-wave resonance condition is satisfied. We verified 



that this is indeed the case by using parametrization (E.l I and the Mathematica package 



Therefore, the main cont ribut ion to the six-wave dynamics arises from the s econd 



(the first non-zero) term in (D.ll 
gives a k-independent result: 



1,2,3 
4,5,6 ' 



This contribution, with the use of (E.l I 



4,5,6 



1,2,3 
4,5,6 



F 2 8/61.4 



(D.2) 



We shall approximate the six- wave interaction coefficient of the IWE| by the constant 



given by (D.2 ). 



E. Locality of the Kinetic Equation Collision Integral 



In this Appendix, we will derive the necessary criteria for the convergence of the colli- 
sion integral in the |KE| We will check for convergence in the limits where a combination 
of wave numbers diverge to infinity or tend to zero. The Dirac delta functions present in 



the collision integral imply that the six- wave resonance condition (A.7) is always satis 



fied. To make our analysis easier we will use a parametrization for the six-wave resonance 
condition. This parametrization is defined by making the values of two wave numbers, 
functions of the remaining four such that the six-wave resonance condition is satisfied. 
This procedure leads to the parametrization: 



k 2 
k 5 



(k 4 -k)(k 3 -k 4 ) 
k + k 3 - k 4 - k 6 
(k 4 -k)(k 3 -k 4 ) 
k + k 3 - k 4 kg 



k + k 3 k 4 . 



(E.la) 
(E.lb) 



We begin by checking for convergence in the IR region of k-space, by determining the 
k-scaling of the |KE| as one wave number, say kg — > 0. Therefore, we must find the |IR| 
k-scaling of each contribution in the collision integral. We will present our analysis for 
a general six- wave |KE| So we assume that the six- wave interaction coefficient has the 
following scaling as fcg — > 0: 

Jim w#; 6 2 cx 4, (E.2) 

kg— >o 



where ^eR. Then in the same limit, the following term in the |KE| behaves as 



lim ?ik n 2 f T'3^4^5 r i6 { n k 1 + 1 + n 3 1 — n 4 1 — n K 1 — n, 1 



k 6 ^o 



oc nkn2n 3 n 4 n5ng (n k 1 + n 2 1 + n 3 1 — n 4 1 



e J (E-3) 
n^ 1 ) oc n 6 oc kg X , 



"5 
—is 



where x is the wave action spectrum exponent for from (66 1. Thus, we can factorize 
out the integral over kg in the |KE| and write it as 



fcg^rig dk 6 oc 2 / x dkg 
Jo 



(E.4) 
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Ultimately, convergence of the |KE| in the |IR| limit corresponds to convergence of integral 
(E.4) as kg — > 0. Hence, convergence of the collision integral is satisfied when 



x < 1 + 2£. 



(E.5) 



Note that all other integrals over k 2 , k 3 , k4, k 5 diverge in exactly the same manner as 
the integral over kg. Moreover, we can check for |IR| divergence if two wave numbers are 
simultaneously small. When two wave numbers on the same side of the six-wave reso- 
nance sextet are small , i.e. k 2 and k 3 or k4 and kg, then we see an integral proportional 



to the square of (E.4) - doubling the rate of divergence. 



If the two wave numbers are on opposite sides of the sextet, then we get convergence 
because of an additional vanishing contribution arising from the term 



lim ( n 

k 2 ,k 6 ^0 



-1 
k 



0. 



(E.6) 



Total convergence and thus locality of the KZ solutions ( 73 ) is only assured if we have 



convergence in both |IR| and UV| limits of k-space. Therefore, we will now proceed in 
calculating the condition of |UV| convergence of the collision integral. 



To check for UV convergence, we will consider the scaling of the collision integral as 
one wave number diverges to infinity. The conservation of momentum, implied by the 
Dirac delta function of wave numbers, means that if one wave number diverges, so must 
a second from the opposite side of the sextelj^] This can be verified by parametrization 
(E.l), i.e. if we force kg to diverge, then so must k 2 . Indeed, due to this second 



divergence, we must integrate over the second diverging wave number and additionally 
compute the k-scaling of the resulting Jacobian. 

To begin, let us assume that the six-wave interaction coefficient scales in the following 
way with the diverging wave number kg: 



lim W^'g oc k%, 



(E.7) 



where rj £ R. Due to the integration of an additional diverging wave number (in this case 
k 2 ), we must consider the k-scaling of the Jacobian corresponding to the transformation 
between kg and k 2 . The Jacobian is expressed as 



fe,2,3 
4,5,6 



dk 2 



2|k + k 3 -k4-k 6 



(E., 



where we have used the fact that k 5 = k + k 2 + k 3 — ki — kg. From equation (E.8), we 



observe that the Jacobian produces a contribution to the collision integral oc k 6 , and 
moreover, we find that the following expression in the collision integral scales as 



lim 

kg— >OQ 



ni tL n 2 n 3 n4n 5 n 6 



oc k 6 (fcg 



(E.9) 



8 Divergence of two wave numbers on the same side of the sextet would violate the Dirac delta function 
involving frequencies. 
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The first term on the right-hand side is a contribution arising from the Jacobian and 
the product of n^s. The second term stems from the difference of n^s (shown in the 
brackets in equation ( 56 1) giving rise to two contributions. The fcg scaling results from the 
difference of the four non-divergent n^s, while the k§ 2 factor comes from the leading order 
Taylor expansion of the two divergent n^s. Therefore, collecting all the contributions 
together, we find that the the condition for the |UV| convergence of the collision integral 
is given by 

r) < x. (E.10) 



Subsequently, the overall locality region for the KZ solutions for the six-wave |KE| is given 
by 

f]<x <l + 2£. (E.ll) 



F. Derivation of the Differential Approximation Model 



In this Appendix we will outline the derivation of the DAM from the KE If one 



assumes that only super-local wave interactions occur, then the |KE| can be simplified into 
a differential model. The |DAM| will describe the evolution of the wave action density 
in w-space - the space determined by angle averaging over wave vectors. This results in 
a dimensional- independent description for n^. Consequently, the wave action density in 
w-space is defined as 

N u = n w - — , 
duj k 

so that the total wave action is equal to 



(F.l) 



duj = / ny dk 



(F.2) 



The KE (561, can then be re-expressed in terms of w. 



dz 



,2,3 



S 4,5,6 n ul n 2 n 3 n ll n 5 n e 
x d^duj^duj^duj^dujQ, 



1 



1 

n 2 



1 

"3 



1 

"4 



1 

"5 



1 

n 6 



where we have defined a new interaction coefficient given by 



S 



w,2,3 
4,5,6 



6 7T 
~6~ 



i 2 A -k.'J.:s\ dk dk 2 dk 3 dk A dk 5 dk 6 

4:5,6 ' du>k duj2 dui 3 duj4 d(jj§ dujg 



(F.3) 



(F.4) 



^4 5 6 has the same symmetry properties as W^'^'g given in (28). The strategy in 



k,2,3 



deriving the DAM from the KE is outlined in 



and involves multiplying (F.3) by 



some arbitrary smooth function, f(uj) = f u , and then integrating with respect to duj. 



Using the symmetry of the interaction coefficient S^^'q, this procedure gives 
n UJ f(iu)u!^ 1 ^ 2 duj = 



S 4,5 2 ,6 3 n LLl n 2 n 3 n 4: n 5 n 6 



x 



1 



11111 

n u n 2 n 3 n 4 n 5 rt 6 

x (fu + h + fa - fi - fa - fa) dui 2 dui 3 dui4duj 5 du! e , 
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(F.5) 



where /, = f(u>i), with i = k, 2, 3, 4, 5, 6. The super-locality assumption of wave interac- 
tions implies that each with i = 2, 3, 4, 5, 6 can be approximated to be close to Uk, such 
that each u>i is within a small deviation of uik, i.e. caj = Wfe(l + Pi) for i = 2,3,4,5,6. 
This permits the Taylor expansion of the two brackets involving n~ s and fs around the 
deviations up to 0(p 3 ). Furthermore, by approximating n^riin-in^a^n^ « n^, and using 
the scale invariance property of S^ 5 ' 6 , we gain the following equation: 



n u f{w)bj 1//2 dui 



So 



dco 2 



duj, 



(F.6) 



where 



So 



P2, 1 +P3, 1+P4, 1 +P5, 1 +Pe) 



X (pl 



P4 P5P6) <*p4d>5,PB 

Then from applying integration by parts we get 



dp 2 dp 3 dp4(ip 5 dp 6 



(F.7) 



hu,f(u>)uj l ' 2 c?w = So 



8uj 2 



^' 2 n 6 , 



doj 2 



f(u) du. 



As /(w) is an arbitrary function, the two integrals must be satisfied for all choices of /(w). 
Therefore, their integrands must equal one another. This results in the formulation of 
the ID AMI 



hu> = S uj 1/2 



duj 2 



u 2 n,. 



d 2 ^ 
du 2 



1 



(F.9) 



G. The Bogoliubov Dispersion Relation 

We will derive the expression for the Bogoliubov dispersion relation for the descrip- 
tion of waves upon a condensate background. The dynamics of waves propagating in 
the presence of a condensate differ from the pure linear waves of the system. The Bo- 
goliubov dispersion relation comprises of a nonlinear correction to the linear frequency 
that emerges due to the existence of a strong condensate. The strategy in deriving the 
formula for the Bogoliubov dispersion relation is to consider the expansion the wave func- 
tion ip(x, z) around a homogeneous condensate background. In the most ideal situation, 
the condensate is represented as the dynamics of the zeroth Fourier mode. However, in 
reality the condensate may be described by a range of Fourier modes situated at low 
k. For simplicity, we shall consider a perturbative expansion of small disturbances on a 
uniform (zeroth mode) condensate. 

The condensate is described by an ^-independent solution to the LWE| Q, i.e. a 
solution of the form ip(x,z) — tp c (z), where tpc(z) is given by 

ip c (z) = ip exp(-iu; c z/2q), (G.l) 

with uj c — ~Io/2ip 2 l 2 and where Iq — \ipo\ 2 is the intensity of the condensate (zeroth 
mode). Solution ip c {z) describes the background rotation of a uniform condensate in Q 
with a rotation frequency of lo c . 
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We expand the LWE is powers of <j>(x, z), the wave function for disturbance upon the 
condensate, and linearize with respect to this perturbation, where we have defined 



TP(x,z)=i> c (z)[l+<t>(x,z)}, (G.2) 



with \<j)(x, z)\ <C 1. Substituting relation (G.2) into the LWE[ and linearizing to 



the first order in (j>(x,z), gives a linear evolution equation for cj){x^z). Assuming the 
disturbance (j>(x, z) takes the form of a single monochromatic plane wave: 

(j)(x, z) = A exp [ ikx — - — — ) + A* exp ( —ikx + - — — ] , (G.3) 
V 2 9 / V 2 ? J 

where A is a complex amplitude of the wave and Qk is the frequency of the waves 
(disturbances) upon the condensate. Then by equating both types of exponentials, we 
can derive the dispersion relation for the propagation of weakly nonlinear waves upon 
the condensate: 

* - J ( 1+ (a4 > 

To obtain the Bogoliubov frequency of the original wave function ip(x, z), we must in- 
clude the frequency in which the condensate is rotating, uj c . Therefore, the Bogoliubov 
dispersion relation for a weakly nonlinear wave in the presence of a condensate is given 
by 




2^H\ V V ^ 2 J ^ 2 l 2 e 



(G.5) 



H. Non-Dimensionalization 



The dimensional models for |lD|[OWT contain several physical constants and param- 
eters. For convenience and clarity, it is preferred to consider the dimensionless forms of 
the equations. Therefore in this Appendix, we present our non-dimensional descriptions 

(|9j. We follow the non-dimensionalization that was 



for the LWE rt7| and the SWE 



performed in [7 1 



ijj = (lp c / y/a)tp* (x*,Z*), X = X c yjc 



z = z c az 



(H.l) 



with tpl = 2K j 'eokgn^l^, z c — 2gZ| and x c = l^. The electrical coherence length of the 
is defined as ^ = v / irK/2Ae{d/V Q ). 



LC 



Non-dimensionalization (H.l), expresses the LWE in the form: 



dip 
oz 



d 2 ^ 
dx 2 



\m 2 



2a ^ dx 2 



(H.2) 



where a is a tunable parameter that adjusts the strength of the second nonlinear term 
and for clarity we have dropped the superscripts on the non-dimensional variables. The 
nonlinearity of the system can be adjusted by increasing or decreasing the magnitude of 
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the wave function tp(x, z), however, parameter a provides additional control of the ratio 
of the the two nonlinear terms. Indeed, the |LWE] was derived in a regime of a long- wave 
limit that implies that (k*) 2 /a<l (where k* is the dimensionless wave number), and 
therefore we must ensure that the long-wave limit it satisfied. 

Similarly, we can derive a dimensionless expression for the SWE using the same 
dimensionless variables given in (H.l|. 
as 

dip 
oz 



This gives the dimensionless version of the SWE 

d- 2 \^\ 2 



d 2 ip 
dx 2 



1 



dx 



(H.3) 



Notice, that in (H.3) a does not appear. 



I. The Intensity Spectrum 



Experimentally, due to the difficulty in measuring the phase of the wave function 
ip(x, z), we cannot determine easily. However, we can measure the spectrum of wave 



Intensity iV k = | 



Theoretically, we are able to relate the fc-scaling for with 



that of the KZ solution derived from the KE (56). In this Appendix, we present the 



derivation of this relationship. First, we must consider the expression for the light inten- 
sity in k-space. We begin by considering the usual definition for the Fourier transform 
in K, this implies 



(H 2 



ip(x)ip* {x)e 



-ika; 



dx. 



(I.la) 
(Lib) 



= J a\a 2 <5 k ' 2 dki dk 2 . 

Hence, for the intensity spectrum N-^ = (|/k| 2 ) this gives 

(|^k| 2 ) = J (a ia2 a* 3 a* 4 ) <5 k ' 4 S^ 2 dki dk 2 dk 3 dk 4 . (1.2) 

The next step is to average over phases, in a |RPA| wave field. This implies that only 
wave number pairings ki — k 4 , k 2 = k 3 and ki = k 3 , k 2 = k 4 will contribute to the 
intensity spectrum. Therefore, the intensity spectrum can be expressed as 



/ (l«i! 2 ) <l«2| 2 ) <5 k <5 k dk! dk 2 + J (K| 2 |a 2 | 2 ) 5\' 2 5\ a dki dk 2 (I.3a) 



mS k dki 



ri\n 2 (5 k ' 2 dk\ dk 2 . 



(I.3b) 



We will examine the intensity spectrum, when the system is in a statistically non- 



equilibrium stationary state and when the KZ solution is realized. Therefore, we can 
assume that the wave action spectrum is of the KZ form, i.e. — Ck~ x , where C is 
a constant determining the amplitude of the spectrum and x is the spectrum exponent. 
The first term of equation (I.3b) contains a Dirac delta function, centered around k = 0. 
This implies that the contribution from this term will only appear at the zeroth mod<r| 



9 Numerically and experimentally this may be seen as a contribution at low wave numbers around 
k = 0. 
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On the other hand, the second term in equation (1. 3b I will contribute on the whole of 



k-space, and will determine the £;-scaling for the intensity spectrum. The fc-scaling of 
the intensity spectrum will only be observed if the integral of the second term is finite. 
We proceed by determining a sufficient condition for its convergence. We can re-express 
the second term by using the Dirac delta function to eliminate one of the integration 
variables, i.e. 



n\ri2 $i' 2 rfki dk2 



7iini_k dki 



= c- 



J fcr x |ki -k|-* dki. 



(1.4) 



If expression (1.4) converges on the KZ solution, then the integral will yield the inten- 



sity spectrum power-law scaling. To check for convergence, we change the integration 
variable to a non-dimensional variable s = ki/k. Then, the intensity spectrum can be 
approximated by the integral: 



C 2 k 



2 h -2x+l I a -*\ a -l\-* ds. 



-2»+l 



-oo 
oo 



x {V 



II 



+ * + l 



ds. 



(1.5) 



Convergence must be checked in the regions where s — > and s —¥ oo. As s — > 0, relation 

(1.6) 



(1.5 1 behaves as 



/ fl -«(| a -l|-* + | s +l| 
Jo 

and therefore, the integral converges for x < 1. 



ds oc / s x ds, 



In the limit when s — > oo, integral (1.5) can be written as 

/•OO 

S - x (\s - l\- x + \s + l\- x ) ds oc / s 



ds. 



(1.7) 



The right-hand side of expression (1.7) converges for x > 1/2. Therefore, the intensity 



spectrum integral, (1.5), is convergent in the region: 

1/2 < x < 1, (1.8) 

where x is exponent of the power-law for the |KZ| solution of the |KE| Experimentally 
at present, we can only produce the inverse cascade scenario and subsequently, we are 
only interested in measuring the intensity spectrum in the inverse cascade regime. The 



exponent of the KZ solution for the inverse cascade is given in equation ( 75b ) [x 



3/5), which lies inside the region of convergence (1.8), (and so the intensity spectrum is 
observable) and corresponds to an intensity spectrum of 



iV k oc k- 1/5 . 



(1.9) 
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